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Damping  and  Trapping  in  2D  Inviscid  Fluids 

N.  Sateesh  Pillai  and  Roy  W.  Gould 

California  Institute  of  Technology,  Pasadena,  California  91125 
(Received  6  July  1993;  revised  manuscript  received  1 1  August  1994) 

We  demonstrate  collisionless  decay  due  to  phase  mixing  of  a  disturbance  in  an  inviscid  2D  fluid 
with  sheared  flow.  Experiments  performed  on  a  cylindrical  pure  electron  plasma,  which  behaves  as 
a  2D  inviscid  fluid,  show  that  a  small  amplitude  quadrupole  excitation  (m  =  2,  ^  =  0  mode)  decays 
exponentially.  At  larger  amplitudes,  the  exponential  decay  is  modulated  by  bounce  motion  of  fluid 
elements  trapped  in  the  azimuthally  traveling  wave  field,  with  the  bounce  frequency  proportional  to  the 
square  root  of  the  excitation  amplitude.  We  also  show  that  the  linear  decay  rate  is  decreased  by  the 
addition  of  external  dissipation.  We  describe  a  calculation  of  the  linearized  response  function. 


PACS  numbers:  47.15.Ki,  52.25.Wz,  52.35.Lv,  52.35.Mw 


Two-dimensional  flow  of  incompressible  inviscid  flu¬ 
ids  has  a  long  history  [1-3].  The  2D  behavior  of  non¬ 
neutral  plasmas,  in  the  drift  approximation,  is  analogous 
to  that  of  inviscid  fluids  [4,5],  with  equipotentials  as  fluid 
streamlines  and  plasma  density  proportional  to  vorticity. 
Very  recently  this  property  has  been  exploited  to  study 
several  classic  problems  in  2D  vortex  dynamics  [5-8]. 
In  sheared  2D  flows,  disturbances  are  convected  in  dif¬ 
ferent  layers  with  different  velocities,  and  this  is  pre¬ 
dicted  to  lead  to  collisionless  decay  of  disturbances  due 
to  phase  mixing.  For  example,  using  an  initial  value  ap¬ 
proach,  Case  [3]  showed  theoretically  that  in  planar  Cou- 
ette  flow  of  an  inviscid  fluid,  small  perturbations  decay 
algebraically  in  the  asymptotic  limit.  Briggs  et  al  [4], 
in  a  theoretical  study,  also  argued  that  small  2D  pertur¬ 
bations  in  a  cylindrical  non-neutral  plasma  column  with 
sheared  flow  are  expected  to  decay  due  to  a  collision¬ 
less  phase  mixing  process  similar  to  Landau  damping  of 
Langmuir  oscillations  [9].  deGrassie  [10]  was  the  first  to 
search  for  this  phenomena  experimentally.  Although  he 
found  the  perturbations  decayed,  the  decay  rates  obtained 
were  about  2  orders  of  magnitude  smaller  than  expected, 
and  they  decreased  with  increasing  amplitude  of  excita¬ 
tion,  i.e.,  a  nonlinear  behavior.  The  sensitivity  of  his 
receiver  was  inadequate  to  observe  the  linear  regime. 

We  present  three  new  experimental  results  on  vortex 
dynamics,  (a)  We  observe  collisionless  damping  of  2D 
disturbances  of  the  form  exp[z>M^  +  ikz  -  i(ot\  with 
m  =  2  and  =  0,  in  an  inviscid  fluid  in  the  linear 
regime,  (b)  We  find  that  as  the  amplitude  of  excitation 
is  increased,  the  decay  rate  is  modulated  periodically 
in  a  manner  consistent  with  bounce  motion  of  trapped 
fluid  elements  in  the  field  of  the  disturbance.  This  is  in 
accordance  with  the  early  picture  of  Kelvin’s  cat’s  eyes 
[2]  and  a  prediction  of  Briggs  et  al  [4]  that  the  bounce 
frequency  should  be  proportional  to  the  square  root  of  the 
perturbed  amplitude.  This  phenomenon  is  qualitatively 
similar  to  the  effect  of  bounce  oscillations  of  particles 
trapped  by  a  Langmuir  mode  [11-13].  (c)  We  show 
that  adding  external  dissipation  to  the  system  reduces  the 


linear  damping  rate,  thus  demonstrating  that  the  mode  has 
negative  energy,  i.e.,  excitation  of  the  mode  reduces  the 
energy  of  the  system. 

We  also  describe  a  calculation  of  the  linear  response 
function  of  a  2D  cylindrical  fluid  which  exhibits  exponen¬ 
tial  decay  at  early  times,  associated  with  a  simple  pole, 
followed  by  a  transition  to  algebraic  decay,  similar  to  the 
asymptotic  prediction  of  Case  [3].  This  response  function 
also  shows  that  the  mode  has  negative  energy. 

Our  experiments  were  performed  using  a  cylindri¬ 
cal  pure  electron  plasma  column  illustrated  in  Fig.  1(a). 
Electrons  emitted  from  a  hot  thoriated  tungsten  filament 
are  injected  into  the  trap  region  where  they  are  con¬ 
fined  in  the  radial  direction  by  a  50  G  magnetic  field 
and  in  the  axial  direction  by  electrostatic  traps  held  at 
-100  V.  The  trapped  electrons  will  rotate  around  the 
cylinder  axis  under  steady  state.  A  radial  monotonically 
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FIG.  1.  (a)  Schematic  of  cylindrical  structure  for  plasma 

trapping  and  excitation,  (b)  phasing  of  first  octupole  for 
exciting  an  m  =  2  disturbance,  (c)  configuration  of  second 
ocutpole  for  signal  reception,  and  (d)  configuration  of  second 
octupole  for  negative  energy  test. 
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decreasing  particle  density  profile  (vorticity  profile  in 
fluid  terms)  produces  a  monotonically  decreasing  angular 
velocity  of  fluid  rotation,  i.e.,  sheared  flow.  The  plasma 
radius  and  lengths  are  1.2  and  40  cm,  respectively.  The 
cyclotron  frequency,  central  plasma,  and  rotation  fre¬ 
quencies  are  about  140  MHz,  10  MHz,  and  350  kHz,  re¬ 
spectively.  The  Reynolds  number  is  estimated  to  be  10"^. 
Other  aspects  of  this  device  have  been  described  earlier 
[14].  The  experiment  is  performed  in  a  repetitive  inject- 
observe-dump  cycle.  The  trapped  plasma  is  excited  by 
a  5  /Lisec  burst  of  300  to  600  kHz  sinusoid  on  one  oc- 
tupole  set.  The  frequency  of  the  rf  pulse  is  adjusted  so  as 
to  maximize  the  response  of  the  w  =  2  mode,  typically 
2-4  times  the  frequency  of  the  well  studied  m  =  1  dio- 
cotron  mode  [9].  As  shown  in  Fig.  1(b),  signal  phases 
on  the  octupole  sectors  generate  an  w  =  2  field.  The 
charge  induced  on  one  sector  of  a  second  octupole  sec¬ 
tion  is  amplified  by  a  low  noise  charge  amplifier  and 
digitized.  The  response  to  the  short  burst,  after  it  has 
ended,  is  a  decaying  sinusoid  as  shown  in  the  inset  of 
Fig.  2.  Each  digitized  signal  is  processed  by  passing  it 
through  a  software  peak  detector  to  obtain  the  envelope 
of  the  decaying  sinusoid.  The  envelopes  of  20  decays 
(typically)  at  the  same  amplitude  are  averaged  and  the 
process  repeated  for  applied  voltage  between  10  mV  to 
1  V.  The  averaged  envelopes  are  shown  on  a  semiloga- 
rithmic  plot  in  Fig.  2.  For  clarity,  the  traces  in  the  plot 
are  stopped  when  the  signal  becomes  comparable  with 
the  system  noise  level.  The  lowest  four  traces  in  Fig.  2 
evidently  represent  the  linear  regime,  since  the  response 
is  simply  proportional  to  the  applied  voltage.  The  lin¬ 
ear  decay  rate  is  y  =  34.3  krad/sec  and  co/y  =  84.  The 


FIG.  2.  Decay  of  the  envelope  of  m  =  2  signal  at  about 
400  kHz,  excited  by  an  rf  pulse  with  voltages  10,  20,  40, 
60,  80,  120,  160,  200,  250,  300,  400,  500,  600,  700, 
and  1000  mV.  Lower  left  curve  corresponds  to  10  mV  rf 
voltage  and  uppermost  curve  corresponds  to  1000  mV.  Plotted 
on  semi  logarithmic  scale  with  vertical  scale  calibrated  in 
dB.  Time  scale  starts  with  /  =  0  at  the  end  of  the  rf  pulse. 
Inset  shows  received  signal  for  500  mV  applied  voltage. 


decay  time  in  the  linear  regime  is  much  shorter  than 
the  electron-electron  collision  time  of  10  msec.  Hence 
this  decay  is  collisionless.  The  initial  decay  rates  for 
all  traces  during  the  first  15  /xsec  are  approximately  the 
same.  Later  in  time  and  at  higher  amplitudes,  the  decay 
rate  is  periodically  modulated,  with  the  shortest  modu¬ 
lation  period  occurring  at  the  highest  amplitude.  At  the 
lower  amplitudes,  the  modulation  period  becomes  larger 
than  the  decay  time  and  the  nonlinear  modulation  phe¬ 
nomenon  disappears.  At  these  low  amplitudes,  a  fluid 
element  cannot  complete  a  trapped  orbit  before  the  mode 
decays  and,  in  effect,  does  not  know  that  it  is  trapped. 
At  higher  amplitudes  a  fluid  element  is  able  to  complete 
a  trapped  orbit  before  the  mode  decays  and  the  modu¬ 
lation  of  the  wave  amplitude  becomes  noticeable.  The 
smallest  observable  bounce  frequency,  94.3  krad/sec,  is 
larger  than  the  linear  decay  rate.  At  the  largest  ampli¬ 
tudes,  the  late-time  decay  is  much  slower,  and  we  find 
a  nonlinear  behavior  in  which  the  decay  rate  decreases 
with  increasing  amplitude,  similar  to  that  described  by 
deGrassie  [10]. 

In  order  to  determine  the  dependence  of  the  bounce 
period  on  the  mode  potential  we  have  analyzed  the  data 
of  Fig.  2  in  the  following  way.  First,  we  determine  a 
least  squares  exponential  fit  to  the  early  part  of  one  of  the 
traces  to  get  the  average  rate  of  decay.  We  then  subtract 
this  exponential  from  that  trace  and  determine  the  bounce 
period  by  measuring  the  time  between  successive  maxima. 
We  observe  that  the  period  lengthens  as  the  mode  decays, 
as  expected.  We  correct  the  mode  potential  for  this  decay 
using  the  exponential  function  and,  assuming  the  initial 
mode  potential  to  be  proportional  to  the  applied  voltage, 
plot  the  resulting  bounce  frequency  versus  the  corrected 
applied  potential.  Figure  3  shows  the  resulting  bounce 
frequency  /  versus  V  on  a  log-log  plot.  A  least  squares 
fit  of  the  form  /  =  AV"  to  the  data,  shown  by  the  solid 
line,  gives  n  =  0.44  and  A  =  2.77  (/  in  kHz  and  V  in 
mV).  Calculations  for  21  other  data  sets,  similar  to  those 
shown  in  Figs.  2  and  3,  give  values  of  n  between  0.42 
and  0.77  with  a  mean  value  of  0.55  ±  0.10.  This  is  equal, 
within  experimental  error,  to  the  value  of  0.5  expected 
for  oscillations  of  trapped  particles  in  the  mode  [4],  and 
reinforces  our  conclusion  that  the  amplitude  modulation 
of  the  decaying  m  =  2  mode  is  due  to  fluid  tiapping  at 
a  radius  where  the  local  angular  velocity  is  equal  to  the 
mode  velocity. 

We  now  sketch  the  theoretical  calculation  of  the 
linear  response  function  of  the  system  associated  with 
a  single  m  value,  namely  m  =  2.  The  response  of 
the  plasma  to  a  burst  of  sinusoid  of  amplitude  V  is 
induced  current  7  in  the  wall  of  the  electrode  (Fig.  1). 
In  the  linear  regime,  the  Fourier  transform  of  I  is 
related  to  the  Fourier  transform  of  V  through  a  transfer 
admittance  Y{co)  :  I{a))  =  V{(o)Y{(o).  It  is  easy  to  show 
that  Y{(x))  =  ia)27T€QLF\F2{rErl 4>)r=^b^  where  L  is  the 
total  length  of  the  trap,  F]  and  F2  are  geometrical  form 
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FIG.  3.  Bounce  frequency  /  versus  voltage  V  (applied  voltage 
corrected  for  the  decay)  plotted  on  a  log-log  scale.  Filled 
boxes  correspond  to  data.  Solid  line  is  a  least  squares  fit 
/  =  AV”,  where  n  —  0.44.  Together  with  21  other  data  sets, 
n  ==  0.55  ±0.1.  Dashed  line  shows  a  line  where  n  =  0.5  for 
comparison. 


factors  (<1)  pertaining  to  length,  mode  number  m,  and 
phasing  of  the  sectors  of  the  input  and  output  octupole 
sections.  Er/(}>  represents  the  ratio  of  the  radial  electric 
field  to  the  potential  and  should  be  evaluated  at  the  wall 
(r  =  b),  Y{(o)  is  just  the  Fourier  transform  of  the  impulse 
response  of  the  plasma.  The  value  of  Er/<t>  at  the  wall 
can  be  obtained  by  solving  the  linearized  equation  for  the 
potential  [4,5]: 


d^(f)  \  d<f>  , 


me 


dno{r)ldr 
r[w  -  mo)o{r)] 


0=0, 

(1) 


where  no(r)  represents  the  particle  density  and  o^oir) 
the  angular  velocity  of  electrons  of  charge  e  confined 
in  a  magnetic  field  Bq,  For  a  band  of  co's  a  mode 
particle  resonance  will  occur  at  a  radius  within  the 
plasma  where  mcooirs)  =  co,  giving  rise  to  a  logarithmic 
singularity  in  0  at  r  =  r5.  To  solve  Eq.  (1)  requires 
a  knowledge  of  no{r)  from  which  a>o(r)  is  determined. 
Power  series  expressions  have  been  obtained  for  special 
profile  functions,  but  for  a  greater  variety  of  profiles  we 
also  integrate  Eq.  (1)  numerically  for  various  w.  In  order 
to  deal  with  the  singularity  at  r  —  in  the  solution  of 
Eq.  (1),  0)  is  assumed  to  have  a  small  positive  imaginary 
part.  Using  this  technique  we  obtain  the  frequency- 
dependent  real  and  imaginary  parts  of  {rEr/<t>)b  that 
are  shown  in  Fig.  4(a).  The  corresponding  envelope  of 
the  impulse  response  is  obtained  by  taking  the  inverse 
Fourier  transform  of  {rEr/(}>)b  and  is  shown  in  Fig.  4(b). 
(r£r/0)6  can  be  fit  rather  well  by  the  expression  for  a 
simple  pole,  plus  a  small  remainder  function  5R(a>),  which 
accounts  for  the  asymmetry:  Y{a))  —  iwBKw  -  a>2  + 
iyi)  +  9^(^e>)(r2  >  0).  The  pole  lies  in  the  lower  half 


Normalized  frequency  (o 


(b) 


FIG.  4.  (a)  Real  and  imaginary  parts  of  irEr/<t>)r~b  as  a 
function  of  frequency  w  for  the  m  =  2.  Central  angular 
velocity  corresponds  to  =  2,  and  angular  velocity  at  the 
wall  {r  =  b)  corresponds  to  o)  =  0.714  for  m  =  2.  The  dotted 
vertical  line  in  the  middle  corresponds  to  the  theoretical 
prediction  of  the  m  =  2  resonant  frequency  cj  —  0.904  for  this 
profile,  (b)  Semi  logarithmic  plot  of  the  envelope  of  inverse 
Fourier  transform  of  {rEr/<f>)r=b‘  Initial  exponential  decay 
change  to  an  asymptotic  algebraic  decay  with  some  interference 
in  between. 


frequency  plane,  a  location  which  can  only  be  reached 
by  analytic  continuation  of  our  solution  of  Eq.  (1),  and 
is  responsible  for  the  early  exponential  decay  (with 
rate  72)  seen  in  Fig.  4(b).  The  precise  way  in  which 
Re{y(cu)}  goes  to  zero  depends  upon  the  way  in  which 
the  density  no('')  goes  to  zero  at  the  plasma  edge,  and 
this  also  determines  the  form  of  the  algebraic  decay  in 
asymptotic  limit  of  the  impulse  response.  The  Q  of 
the  resonance  W2f2y2  and  the  frequency  0)2^  relative  to 
the  central  rotation  frequency,  depends  on  details  of  the 
density  profile.  Some  profiles  exhibit  a  Q  close  to  the 
experimentally  observed  value  of  84.  The  form  of  y(<y) 
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is  qualitatively  similar  for  the  various  density  profiles  and 
modes  that  we  have  studied  numerically. 

An  important  feature  of  the  admittance  calculation 
is  that  the  constant  B  is  negative,  opposite  from  that 
expected  for  passive  systems  (e.g.,  a  tuned  RLC  circuit). 
Neutrally  stable  (undamped)  negative  energy  modes  are 
destabilized  by  removal  of  energy  through  dissipation. 
This  has  been  well  documented  for  the  /«  =  1  diocotron 
mode  [15].  In  our  case,  the  m  =  2  mode  is  already 
damped,  though  not  by  dissipation  but  one  might  expect 
that  the  addition  of  dissipation  would  reduce  the  damping 
rate  and  perhaps  even  turn  it  into  a  growing  mode. 
Using  the  admittance  function  one  can  show  that  a 
resistive  or  dissipative  boundary  condition  at  the  outer 
edge  of  the  plasma  reduces  the  growth  rate.  When  a 
resistor  R  is  placed  in  parallel  with  the  octupole  section, 
as  in  Fig.  1(d),  the  following  equation  holds:  Via))  + 
(1//?)  =  0.  Using  the  simple  pole  approximation  and 
neglecting  the  small  remainder  function  ?)t(<i>)  we  obtain 
for  mode  frequency  co  —  a)2  -  iiyi  +  coiBR),  The  last 
term  shows  that  when  B  <  0  there  should  be  a  reduction 
in  the  damping  rate  of  the  mode  proportional  to  the 
resistance  /?,  whereas  when  ^  >  0  the  damping  rate  is 
increased. 

External  dissipation  is  easily  added  experimentally  by 
connecting  resistors  to  unused  segments  of  the  octupole 
sections.  Two  precautions  are  necessary.  First,  the  resis¬ 
tors  must  be  connected  so  as  to  affect  the  m  =  2  mode  but 
not  the  m  =  1  mode,  because  the  latter  would  be  imme¬ 
diately  destabilized.  This  is  accomplished  connecting  the 
resistor  simultaneously  to  two  diametrically  opposite  seg¬ 
ments  of  an  octupole,  as  shown  in  Fig.  1(d).  Second,  the 
cable  capacitance  must  be  “tuned  out”  by  adding  a  paral¬ 
lel  inductor,  selected  so  that  parallel  resonant  frequency  is 
equal  to  the  frequency  of  the  m  =  2  mode.  Experimen¬ 
tally  it  is  found  that  the  connection  of  an  external  resistor 
to  the  wall  electrodes  in  this  fashion  causes  the  damping 
rate  to  decrease  rather  than  increase  and  with  high  enough 
resistance;  the  mode  is  seen  to  be  destabilized.  The  re¬ 
duction  in  decay  rate  by  external  dissipation  is  shown  in 
Fig.  5.  Theoretically,  since  the  form  of  Y{(i))  is  similar 
for  higher  modes  also,  one  expects  similar  behavior  for 
the  m  =  3  and  higher  modes. 

In  summary,  our  experimental  observations  reveal,  for 
the  first  time,  three  important  characteristics  of  an  aw  =  2 
fluid  perturbation.  First  is  collisionless  damping  of  small 
perturbations  with  evidence  for  a  dominant  pole  which 
follows  from  the  observed  exponential  decay  of  the  linear 


FIG.  5.  Decay  rate  of  the  m  =  2  response  versus  resistance 
when  a  resistor  is  connected  to  the  wall  electrode.  Line  is  a 
least  squares  fit  of  the  experimental  data  obtained  from  a  large 
number  of  responses. 

response.  Second  is  evidence  at  large  amplitudes  for 
fluid  trapping  by  the  mode,  with  a  bounce  frequency 
proportional  to  the  square  root  of  the  amplitude.  Third, 
the  decaying  disturbance  has  negative  energy  and  decays 
less  rapidly  when  external  dissipation  is  added. 

This  work  was  supported  by  ONR  Grant 
No.  N00014-89-J-1264. 
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This  paper  focuses  mainly  on  the  dynamics  of  two-dimensional  cylindrical 
nonneutral  plasma.  After  reviewing  some  highlights  of  the  nonneutral  plasma 
dynamics,  some  recent  2D  results  are  described:  vortex  dynamics,  diocotron 
instabilities  of  hollow  profiles,  collisionless  damping  of  modes  and  fluid 
trapping  by  modes,  fluid  echos,  the  cyclotron  center  of  mass  modes  and  warm 
plasma  Bernstein  modes,  their  use  for  temperature  determination.  We  call 
attention  to  some  unsolved  problems. 
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I.  INTRODUCTION. 


A  brief  chronology  of  developments  in  nonneutral  plasmas  is  shown  in 
Table  I.  Work  on  nonneutral  plasmas  goes  back  at  least  to  the  pioneering 
paper  by  Leon  Brillouin,*  in  which  he  first  described  steady  state  flows  of 
charged  particle  beams  in  a  magnetic  field,  a  form  of  magnetic  focussing  in 
which  the  repulsive  space  charge  electric  field  E  generated  by  charged 
particles  is  balanced  by  a  Lorentz  v  x  B  force  of  the  particles  moving  through 
a  static  magnetic  field.  He  described  both  planar  and  cylindrical  steady  state 
flows,  illustrated  in  Fig.  1  for  the  case  of  electrons.  Both  were  later  used 
for  magnetic  focusing  of  electron  beams  in  various  types  of  microwave  devices. 
One  of  the  earliest  attempts  to  understand  the  dynamics  of  such  beams  was  by 
McFarlane  and  Haye  .  In  attempting  to  explain  the  operation  the  highly 
successful  magnetron  device,  they  studied  the  behavior  of  small  perturbations 
in  a  planar  Brillouin  beam  and  showed  that  they  could  be  unstable.  Because  of 
the  sheared  nature  of  the  flow,  the  instability  was  termed  the  slipping  stream 
instability  and  also  the  diocotron  instability.  However,  attempts  to  give  a 
quantitative  theory  of  the  magnetron  were  not  very  successful,  but  Brillouin 
flow  found  important  applications  in  beam  type  microwave  devices,  such  as  the 
traveling  wave  tube^  (cylindrical  beams)  and  the  French  M-type  backward  wave 
oscillator^  (planar  beams).  In  these  applications  good  quantitative  theories 
of  the  of  the  beam  dynamics,  together  with  their  interaction  with  slow-wave 
circuits,  were  possible.^  In  the  case  of  the  M-type  backward  oscillator, 
the  instability  of  the  thin  beam  resulted  in  a  reduced  start-oscillation 
current.  In  the  case  of  magnetically  focused  hollow  beams,  the  instability 
resulted  in  beam  breakup.  However,  most  microwave  device  researchers  of  the 
time  did  not  think  of  their  beams  as  plasmas. 


The  science  of  plasma  physics  began  to  prosper  in  the  1960s  because  of 
the  interest  in  fusion  energy,  accelerators,  and  ion  sources.  Work  on 
nonneutral  plasmas  accelerated  in  the  1970s.  The  cylindrical  Brillouin 
equilibrium  was  verified  over  the  entire  range  of  rotational  velocities.® 
Charged  particles  were  injected  and  trapped  in  magnetic  mirror  fields.®  A 
variety  of  configurations  were  devised  and  studied,  both  experimentally  and 
theoretically.  Considerable  attention  was  given  to  the  equilibria  and  their 
stabi li ty . ^ ^  The  diocotron  instability  was  further  analyzed  and  the 
analogy  between  with  2D  incompressible  inviscid  fluid  flow  was 
recognized. An  excellent  and  comprehensive  review  of  this  phase  of  the 
nonneutral  plasma  research  is  found  in  the  two  volumes  on  nonneutral  plasmas 
by  Davidson. 

Confinement  of  pure  electron  and  pure  ion  plasmas  using  Penning  traps 
became  important  in  the  1970s  and  developed  along  two  tracks,  largely 
independent  at  first.  Both  configurations  employ  electrostatic  fields  for 
axial  confinement  and  are  illustrated  in  Fig.  2. 

The  Penning  trap  with  hyperbolic  electrodes^”^  was  introduced  to  contain 

charged  particles  for  precision  spectroscopic  measurements.  In  hyperbolic 

traps  the  electrodes  are  shaped  so  that  the  confining  potential  is 
2  2 

proportional  to  r  -2z  .  This  makes  the  orbit  of  a  single  particle  in  the  trap 
particularly  simple.^®  In  spectroscopic  work,  plasma  effects  are  undesireable 
and  minimized.  There  has  been  very  extensive  work  by  the  chemists  in 
ion-cyclotron  resonance  mass  spectrometry  in  traps.  Positrons  and  antiprotons 
in  sizeable  numbers  are  trapped,  contained,  and  transported. The 
densities  in  hyperbolic  traps  can  be  made  sufficiently  high  that  plasma 
effects  become  important  and  very  low  temperatures  have  been  acheived  by  laser 


cooling.^®  At  low  temperatures  the  thermal  equlibrium  in  a  quadratic  trap 

potential  is  a  rigidly  rotating  spheroid  of  constant  density.  These  states 

are  routinely  observed.  Dubin^®  has  given  a  theory  of  the  plasma  modes  of 

cold  spheroidal  plasma  in  the  rigid  rotor  state.  These  modes  are  now  used  as 

20  21 

a  diagnostic  to  determine  plasma  densities,  shapes,  and  temperatures,  ’ 
the  latter  with  an  extension  of  Dubin’s  theory. 

Malmberg  and  coworkers^^  developed  the  cylindrical  Penning  trap  to  to 
study  basic  plasma  phenomena  in  a  relatively  idealized  situations.  Many 
phenomena  have  been  identified  and  studied.  >  >  »  »  One  of  the  major 

developments  to  come  from  this  line  of  investigation  is  the  new  experimental 
method  for  the  study  of  2D  vortex  dynamics,  a  subject  of  longstanding 
interest  in  the  field  of  fluid  mechanics.  A  good  summary  of  more  recent 

research  on  nonneutral  plasmas  will  be  found  in  two  A.I.P.  Conference 

„  ,.  28,29 

Proceedings. 

The  remainder  of  this  paper  deals  with  2D  phenomena  in  cylindrical  pure 

electron  plasmas  and  is  organized  as  follows.  In  Sec.  II  we  summarize  the 

Brillouin  equilibrium,  the  2D  drift-dynamic  equations  for  low  frequencies,  low 

densities,  and  zero  temperature.  We  review  their  similarity  with  the  Euler 

101123 

equations  for  2D  inviscid  imcompressible  fluid  flow  ’  ’  and  cite  a  few 

simple  but  important  results.  In  Sec.  Ill  we  discuss  the  linearized  equations 
for  situations  in  which  the  density  and  angular  velocity  depend  on  radius 
(shear-flow)  and  instabilities  when  the  density  and  angular  velocity  profiles 
are  nonmonotonic.  In  Sec.  IV  we  discuss  internal  perturbations,  collisionless 
(Landau-like)  damping  of  the  response  to  an  applied  pulse.  These  can  be 
understood  in  terms  of  phase  mixing  because  of  the  different  angular  velocites 
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with  which  persistent  density  perturbations  are  convected  at  different  radii. 
Experimental  evidence  is  presented  for  trapping  of  the  fluid  by  the  wave.  In 
Sec.  V  we  discuss  the  possibility  of  echos  which  arise  out  of  the  nonlinear 
interaction  between  two  applied  pulses.  In  Sec.  VI  we  discuss  high  frequency 
(cyclotron)  modes  where  inertial  terms  and  plasma  temperature  are  important. 
Finally  in  Sec.  VII  we  enumerate  some  unsolved  problems  connected  with  the 
above . 

II.  2D  CYLINDRICAL  DYNAMICS  AND  THE  INVISCID  FLUID  ANALOGY. 

The  Brillouin  steady  state  for  a  cold  cylindrical  nonneutral  plasma  is  one 
in  which  the  inward  magnetic  force  due  to  the  rotational  angular  velocity 
cancels  the  outward  electrostatic  repulsive  force  of  the  charges  and  the 
outward  centrifugal  force  at  every  radius.  When  the  density  n^(r)  is 
independent  of  radius,  the  angular  velocity  of  rotation  is  also 

independent  of  radius,  the  so-called  rigid  rotor  state.  The  rotation 
frequency  is  determined  by  the  force-balance  equation 

(jj  ^  -  (jj  w  +0;  ^/2  =  0  (1) 

0  0  c  p 

and  the  solutions  of  this  equation,  which  are  well  known,  are  depicted  in 

8  2 

Fig.  3,  along  with  experimental  data.  At  the  maximum  density,  = 

2 

/2,  the  Brillouin  limit.  Most  experiments  on  pure  electron  plasma,  which 
will  be  the  focus  of  this  paper,  are  at  low  density,  in  the  lower  left  corner 

2  2 

of  Fig.  3,  where  .  Small  amplitude  2D  perturbations  from  the  rigid 

rotor  steady  state  are  surface  waves,  traveling  around  the  perimeter  ~ 
exp( imtf-iwt) ,  as  illustrated  in  Fig.  4.  There  are  discrete  low  frequency  and 
high  frequency  surface  modes: 

a;  =  [m  -  1  +  (a/b)^"*]  w 

0 


(2a) 
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w  =  w  +[in-l  -  (a/b)^™]  w  (2b) 

c  0 

where  a  and  b  are  the  radii  of  the  plasma  and  a  surrounding  metal  cylinder, 
respectively.  The  m  =  1  modes  are  simple  center  of  mass  modes  and  the  sum  of 
their  frequencies  is  The  frequency  of  the  m  =  1  low  frequency  mode,  w  = 

w^Ca/b)*,  depends  only  on  the  mean  density  within  the  cylinder.  These  modes 

have  no  z-dependence  (k^  =  0).  With  z-dependence  ~  exp(ikz)  the  modes 

propagate  in  the  z  direction  as  well.  Axially  propagating  modes  have  been 

analyzed  ’  ’  and  they  are  similar  to  the  Trivelpiece-Gould  modes  of 

a  neutral  plasma,  except  that  plasma  rotation  adds  a  new  feature.  For  m  =  0, 

plasma  rotation  does  not  matter,  and  the  axially  propagating  modes  have  been 

2  2 

studied  experimentally  . 

For  low  density  and  low  frequencies  the  inertial  terms  can  be  neglected, 
e.g.  the  first  term  in  Eq.  (1),  and  the  drift  approximation  can  be  used,  and 
the  electron  fluid  velocity  is  given  by  v  =  (E  X  B  )/B„.  For  electrostatic 
disturbances  E  =  -  so  that  2D  fluid  flow  is  incompressible,  V*v  =  0.  Since 

^  =  -V4>  X  B^/B^,  the  flow  is  along  equipotentials  and  <f>  is,  aside  from  a 
factor  1/Bp  the  stream  function.  Furthermore  the  vorticity  ( z-component)  V  X 

_v  =  V  ^  is,  for  a  pure  electron  plasma,  equal  to  ne/e^.  The  density  is 
related  to  the  velocity  by  the  continuity  equation  which,  because  V*_v  =  0,  can 
be  written  as  dn/3t+v«Vn  =0.  This  means  that  the  density  is  simply  convected 
with  the  flow.  The  fact  that  these  equations  have  the  same  form  as  the  Euler 
equations  for  an  inviscid  ordinary  (uncharged)  fluid  was  recognized  early  in 
the  study  of  non-neutral  plasmas. A  comparison  of  the  the  two  sets  of 
equations  is  given  in  Refs.  11  and  24.  However  it  is  only  in  the  past 
half-dozen  years  that  this  feature  has  been  exploited  experimentally  and  it 


has  produced  some  beautiful  experimental  results  in  2D  vortex  dynamics , 
and  made  contact  with  a  large  body  of  literature  in  theoretical  and 
computational  vortex  dynamics  in  the  field  of  fluid  mechanics. 

Perhaps  one  of  the  most  striking  examples  is  the  merger  of  two  vortex 
patches  when  their  initial  separation  is  less  than  a  critical  value,  about  1.7 
times  the  diameter  of  the  vortex  patch.  This  was  a  long  standing  prediction 
in  fluid  mechanics  but  it  was  first  demonstrated  experimentally  in  a 
definitive  manner  by  Fine  et  al  in  a  pure  electron  plasma.^®  This  work  is 
also  an  excellent  example  of  the  ability  to  image  the  line-averaged  plasma 
density  (proportional  to  vorticity)  when  it  is  dumped. 

III.  RADIAL  DENSITY  PROFILES  AND  STABILITY. 

The  instability  of  the  planar  beam  described  in  the  introduction  has  a 
cylindrical  analog.  It  was  first  demonstrated  experimentally  in  hollow 
electron  beams  by  Kyle  and  Webster,^  when  they  showed  that  as  the  hollow 
beam  propagates  it  tends  to  break  up  into  a  number  of  spiral  vortices,  the 
number  depending  on  the  thickness  of  the  beam.  This  is  the  simplest  example 
of  a  non-monotonic  density  profile:  constant  density  between  two  radii,  r^  and 
Tj*  zero  for  r  <  r^  and  for  r  >  Tj-  incompressible  and 

sheared,  with  the  angular  velocity  of  electrons  at  the  outer  surface  greater 
than  at  the  inner  surface.  The  stability  boundaries  for  this  geometry  were 
determined  by  Levy.^®  The  notion  of  surface  waves,  running  around  the  inner 
and  outer  surfaces  at  different  velocities  is  useful.  The  wave  on  the  outer 
surface  has  negative  energy  and  the  wave  on  the  inner  surface  has  positive 
energy.  If  the  beam  is  not  to  thick  (compared  with  the  azimuthal  wavelength) 
the  fields  of  the  positive  and  negative  energy  waves  overlap  and  their 


velocities  are  suffiently  close  that  they  are  coupled  and  this  allows  them 
both  to  grow  exponentially.  It  is  easy  to  show  that  the  growth  rate  depends 
on  the  ratio  of  the  beam  thickness  to  azimuthal  wavelength  and  has  a 

o  o 

maximum,  ’  as  shown  in  Fig.  5.  In  experiments  where  the  instability  is 
initiated  by  noise,  the  wavelength  with  the  fastest  growth  rate  should  be 
observed  and  this  decreases  with  beam  thickness.  This  feature  of  instability 
has  been  nicely  demonstrated  by  Peurrung  and  Fajans  in  a  cylindrical  trap 
as  shown  in  Fig.  6. 

In  most  pure  electron  plasmas,  the  radial  density  distribution  is  not 
uniform,  though  approximately  axi-symmetric .  In  this  case  the  angular 
velocity  also  depends  on  radius,  i.e.  the  flow  is  sheared.  This  significantly 
increases  the  complexity  and  subtlety  of  the  dynamics.  The  angular  velocity 
becomes  a  function  of  radius  and  is  determined  by  the  density  through  the 
equation 

-r 

u  (r)  = — - — -  n„(r')r'dr'  (3) 

0  0 

It  is  convenient  to  normalize  the  density  and  angular  velocity  profiles  to 
their  central  values  s(r)  where 

u  (0)  =  u^i0)/2u  is  the  central  rotation  frequency.  Such  profiles  are 
0  p  c 

expected  to  evolve  to  the  rigid  rotor  steady  state  but  the  time  to  do  so  in 
many  experiments  is  too  long  for  this  to  be  achieved.  This  evolution  is 
sufficiently  slow  that,  for  most  oscillation  and  wave  phenomena,  the  profiles 
can  be  regarded  as  steady.  Examples  of  monotonic  and  nonmonotonic  density 
profiles  and  the  resulting  angular  velocity  profiles  which  will  be  discussed 
in  this  paper  are  shown  in  Fig.  7.  The  last  of  these.  Fig.  7d,  is  the  hollow 
beam  described  briefly  above. 


For  this  paper  we  are  interested  in  the  smooth  profiles  shown  in  Fig.  7b 
and  Fig.  7c.  Linear  inodes  of  these  profiles  have  many  similarities  (and  some 
differences)  with  longitudinal  plasma  waves  in  a  Vlasov  plasma.  In  Ref.  11 
smooth  profiles  were  examined  and  it  was  shown  that  a  necessary  condition  for 
instability  (existence  of  a  temporally  growing  eigenmode)  of  such  a  system  is 
that  dn^(r)/dr  change  sign,  i.e.  that  the  profile  is  nonmonotic,  such  as  in 
Fig.  7c.  This  is  not  a  sufficient  condition  however  and  specific  profiles 
have  to  be  examined.  Fig.  7c  depicts  (approximately)  the  profile  used  in 
the  experiment  of  Driscoll  et  al  which  demonstrated  that  the  m=2  mode  was 
unstable.  In  that  work,  a  small  seeded  m=2  perturbation  was  seen  to  grow 
exponentially  from  the  linear  regime  and  to  form  two  vortices  (density  clumps) 

O  A 

which  eventually  merge,  as  illustrated  in  Fig.  8.  Eventually,  the  high 
density  annulus  ends  up  in  the  center  and  the  density  profile  is,  in  a  coarse 
grain  sense  mononotonical ly  decreasing. 

3  4 

This  situation  was  also  studied  with  a  2D  fluid  code  which  solved  the 
nonlinear  equations  described  in  Section  II  on  a  256  x  256  square  grid  on 
which  the  potential  at  all  points  outside  an  inscribed  circle,  corresponding 
to  the  conducting  wall,  was  set  to  zero.  A  successive  overrelaxation  Poisson 
solver  was  used  and  the  density  was  advanced  in  time  using  a  15  point  Arakawa 
algorithm.  The  density  at  6  times  in  the  evolution,  corresponding  to  about  20 
rotation  times  of  the  column  is  shown  in  Fig.  9.  The  first  "frame",  where  the 
m=2  perturbation  of  the  annulus  is  hardly  visible,  is  in  the  linear  regime. 

By  the  third  frame  the  two  large  vortices  have  formed  and  they  are  beginning 
to  merge  in  the  fourth  frame.  Two  density  holes  persist  near  the  edge  for 
some  time.  Recently,  Huang  has  shown  that  the  asymptotic  experimental  density 
profile  is  reasonable  well  fit  by  a  minimum  enstrophy  theory  rather  than  a 
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maximum  entropy  theory. 
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IV.  CONTINUUM  OF  MODES.  OOLLISIONLESS  DAMPING.  AND  FLUID  TRAPPING. 

Before  describing  experimental  results  on  collisionless  damping  and 
trapping  we  review  the  linear  behavior  of  disturbances  in  stable  profiles, 
following  an  approach  similar  to  that  Ref.  11.  We  use  the  usual  notation  for 
unperturbed  and  perturbed  quanties,  denoting  them  by  subscripts  0  and  1, 
respectively,  and  assume  angular  dependence  of  exp(imd).  The  linearized 
Poisson’s  equation  becomes 
.2 


d>.  ,  d(^,  2 

=  „  e/E 

iv^  'dr  r'  ‘  ‘  " 


(4) 


and  the  linearized  continuity  equation  becomes 


an, 

at 


+  imw  n  - 


dn. 
im  0 


o“i  rB„  dr  '>'1 
0 


=  0 


(5) 


r*" 

(j)  (r)  =  -e  G  (r,r')n  (r')27rr'dr'  +  <j> 

1  j  _  m  1  1 


ext 


(6) 


2 

where  we  have  used  the  drift  equation  X  /  B^,  the  fact  that 

V  =  rw  is  in  the  6  direction,  and  neglected  the  nonlinear  term  v  *7  n, . 

0  0  7  0  J  I 

n^(r)  is  assumed  to  be  known  and  ^^^(r)  is  determined  from  it  by  Eq.  (3).  One 
can  write  the  solution  of  Eq.  4.  as 

'o 

where  G  (r,r')  is  the  Green’s  function  for  the  problem  with  ^,(0)  =  <j>Ah)  =  0. 
m  11 

In  Eq.  (6)  the  potential  is  split  into  the  part  produced  by  the  plasma  and  the 
externally  applied  part.  4>  .  is  the  potential  due  to  the  voltages  applied 

1  6X  u 

to  the  conducting  wall. 

There  are  two  types  of  problems  which  one  might  want  to  solve;  (a)  an 
initial  value  problem  where  the  initial  density  n^(r,0)  is  given  and  one  wants 
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to  know  how  the  density,  potential,  etc  evolve  subsequently  in  time  or,  (b)  a 
situation  where  a  voltage  pulse  is  applied  to  a  segment  of  the  conducting 
electrode  at  t=  0  when  the  plasma  is  otherwise  undisturbed,  and  one  wants  to 
know  how  the  plasma  responds  to  the  pulse.  A  Laplace  transform  with  respect 
to  time  is  suitable  for  either  type  of  problem.  In  the  latter  case,  n^(r,0)  = 
0  and  Eqs.  (4)  and  (5)  become 

d^$,  ,  d$,  2 


dn 

= -TT  dr*i 

0 

where  $^(r,a;)  and  Nj(r,u;)  are  the  Laplace  transforms  of  the  potential  and 
density.  We  use  for  the  Laplace  variable  s  =  icj  with  u  in  the  upper  half 
plane.  Eliminating  gives 

d^$  d$  2  dn^/dr 

- L  +  i__L_I!L$  +  _e]n _ 2 _  .  ^  , 

i_2  r  ,  2  ^  rB  €  (w-mw  ) 

dr  dr  r  o  o  o 

This  equation  has  been  extensively  studied.  Briggs  et  al^^  were  the  first 

to  give  a  general  discussion  of  the  analytic  properties  of  in  the  complex  w 

plane  and  to  note  the  similarity  of  this  situation  with  that  of  Landau  damping 

of  longitudinal  plasma  waves.  $^(r,w)  and  N^(r,w)  are  defined  in  the  upper 

half  u  plane  and  as  the  real  w  axis  is  approached  from  above,  the  denominator 

of  Ekj.  (9)  becomes  singular  if  the  real  part  of  w  lies  in  the  range 

m  .  <  Re{w}  <  m  w  ,  (10) 

omin  omax  ’ 

where  w  .  and  are  the  minimum  and  maximum  rotation  frequencies  in  the 

omin  Omax  ^ 

plasma.  This  gives  rise  to  a  branch  line  along  the  real  frequency  axis 
between  these  limits  as  shown  in  Fig.  10a.  The  potential  function  can  be 


analytically  continued  into  the  lower  half  plane  as  discussed  in  Ref.  11, 
leading  to  the  deformed  contour  of  Fig.  10b. 

As  with  Langmuir  waves,  one  can  show  that  Eq.  (9)  has  a  continuum  of 
singular  eigenmodes  (generalized  functions)  in  the  frequency  range,  similar 
to  the  Van  Kampen-Case  modes  ’  for  Langmuir  oscillations.  This  leads  to 
two  different  views,  or  descriptions,  of  collisionless  damping  which  are 
equally  valid.  Although  the  diocotron  mode  problem  is  similar  in  many  ways  to 
the  Landau  problem,  it  is  different  in  two  significant  ways:  the  frequency 
range  is  limited  by  Eq.  (10)  and  each  (angular)velocity  is  associated  with  a 
different  radius. 

3  6 

For  comparison  with  experiment  we  are  primarily  interested  in  the 
response  of  the  plasma  to  an  applied  pulse.  We  need  solutions  of  Eq.  (9)  in 
which  $^(b,a;),  the  Laplace  transform  of  the  potential  at  the  conducting  wall 

is  the  applied  potential.  must  vary  as  r™  near  the  origon.  To  compare  with 
experiment  we  need  the  radial  electric  field  at  r=b  from  which  the  current 
flowing  to  the  wall  can  be  obtained.  For  a  parabolic  density,  f(r)=  1  - 
(r/a)^,  Corngold^®  has  shown  that  the  solution  of  Eq.  (9)  can  be  expressed  in 
terms  of  the  hypergeometric  function  and  the  location  of  the  pole  illustrated 
in  Fig.  10b  can  be  determined  analytically. 

We  are  interested  in  more  general  profiles,  n^(r)  and  ^^(r),  where 
solutions  are  not  necessarily  expressible  in  terms  of  known  functions.  Our 
approach  is  to  numerically  integrate  Eq.  (9)  for  a  specific  density  profile 
and  for  a  discrete  set  of  equally  spaced  frequencies  along  the  contour 
shown  slightly  above  the  real  frequency  axis  in  Fig.  10a.  We  obtain  the 
temporal  response  (wall  current)  from  an  inverse  Discrete  Fourier  Transform 


(DFT)of  I ( r /^r  I  .  To  compare  with  experimental  measurements  the 
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admittance  function  is  defined  by 

Y(w)  =  +iu;27rf^F^F2|r[d$^(r,a;)/dr]/$^(r,u;)  ,  (11) 

where  F^  and  F^  are  sector  form  factors.  The  admittance  function  is  the 
Fourier  transform  of  the  response  to  current  in  sector  2  due  to  an  impulse 
applied  to  sector  1.  For  smoothness  of  this  function,  there  is  a  requirement 
that  the  ratio  of  the  step  size  in  the  radial  integration  and  the  step  size  in 
frequency  satisfy  the  inequality  Ar  «  Aw/(dw^/dr).  Evaluation  of  the 
function  Y(w)  at  frequencies  at  a  distance  tr  slightly  above  the  real  frequency 
axis,  followed  by  taking  its  inverse  with  a  DFT  which  assumes  u  to  be  real, 
gives  rise  to  an  additional  exponential  damping  factor,  which  can  easily 

be  corrected  for.  The  real  and  imaginary  parts  of  the  the  principal  factor  in 
the  admittance  function,  /dr ] are  shown  for  a  typical  monotonic 

density  profile  and  for  m=2  in  Fig.  11.  Also  shown  in  Fig.  11  is  the  m  =  2 
impulse  response  obtained  from  the  inverse  DFT  of  YCw) .  The  admittance  has  a 
nonzero  real  part  only  in  the  the  frequency  range  described  by  Eq.  (10)  when 
there  is  a  resonance  in  the  denominator  of  Eq.  (9).  The  admittance  function 
is  defined  so  that  Re{Y}  >  0  corresponds  to  absorbtion.  Since  we  find  that 
Re{Y}  <0  in  this  range,  energy  is  released  when  a  disturbance  is  excited,  so 
the  system  is  a  negative  energy  system. 


The  real  and  imaginary  parts  of  rE^/$  shown  in  Fig.  11  can  be 
approximated  by  a  simple  pole  below  the  real  frequency  axis  plus  a  small 
remainder  function 


Y(w) 


A 

-i(w-a;2+i72) 


+  ??(a>) 


(12) 


where  A  is  negative  because  of  the  negative  energy  property,  Eq.  (12)  is 
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consistent  with  Fig.  10b,  which  shows  a  simple  pole  below  the  real  freqency 
axis.  R(a;)  comes  from  integrals  along  the  branch  cuts.  The  exponential  decay 
which  dominates  the  early  time  behavior  is  associated  with  the  pole  and  the 
algebraic  decay  which  appears  later  is  associated  with  5?(a;) .  While  the 
precise  details  of  Y(a;),  and  thus  of  the  impulse  response,  are  profile 
dependent,  the  general  features  are  the  same  for  qualitatively  similar 
profiles.  The  frequency  w  corresponds  to  twice(m=2)the  rotational  angular 
velocity  at  a  radius  near  the  outer  limit  of  the  plasma  (~0.8a). 

Collisionless  damping  of  the  m  =  2  mode  in  the  linear  regime  has  recently 
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been  observed  by  Pillai  using  the  experimental  apparatus  depicted  in 
Fig.  12.  The  envelope  of  the  decaying  sine  wave,  which  is  excited  by  a 
5/isec  sinusoidal  burst  (about  3  cycles),  is  shown  in  Fig.  13  for  various 
amplitudes  of  the  burst.  Each  trace  is  the  average  of  20  "identical"  traces 
and  an  example  of  the  actual  sinusoidal  signal  is  shown  in  the  inset.  The 
lowest  three  traces  are  essentially  in  the  linear  regime.  As  the  amplitude  of 
the  excitation  is  increased,  "bounce"  oscillation  begin  to  appear,  indicating 
trapping  of  the  fluid  near  the  resonant  layer.  An  analysis  of  the  bounce 
frequency,  when  corrected  for  decay  of  the  wave  amplitude,  shows  that  it  is 
proportional  to  the  square  root  of  the  applied  voltage.  Trapping  in  fluids  in 
shear  flow  was  anticipated  by  Kelvin  with  his  "cats-eyes"  describing  fluid 
trajectories.  In  Ref.  11  it  was  predicted  the  "bounce"  or  trapping  frequency 
should  be  given  by 

w,  «  m  V|$,w'(r)/r|  (12) 

b  '10  'r^ 

o  £ 

Pillai  observes  that  the  bounce  frequency  is  proportional  to  the  applied 
potential  to  the  power  0.55  ±  0.10  using  22  data  sets.  Using  the  2D  fluid 


code  described  earlier,  both  collisionless  damping  and  amplitude  modulation  of 
the  decaying  wave  due  to  trapped  particles  have  been  verified^^  Fig.  14. 
shows  a  calculation  using  this  fluid  code  of  the  charge  induced  on  one  octant 
of  the  conducting  cylinder  when  a  three  cycle  sinusoidal  burst  is  applied. 

The  linear  decay  rate  agrees  well  with  that  obtained  by  numerical  solution  of 
the  linear  differential  equation  and  the  bounce  frequency  is  proportional  to 
the  square  root  of  the  amplitude.  deGrassie  first  looked  for  linear 
collisionless  damping,  but  his  detection  sensitivity  did  not  allow  him  to 
study  the  linear  regime.  Instead  he  found  a  much  slower  decay  whose  rate 
decreased  with  increasing  amplitude,  much  like  the  late  time  behavior  in  Fig. 
13.  He  also  failed  to  see  the  bounce  motion. 

The  negative  energy  feature  of  the  perturbations  has  an  interesting 
consequence.  When  energy  is  removed  from  the  system,  by  external  resistors 
for  example,  the  amplitude  should  grow.  This  was  previously  established  for 
the  m=l  center  of  mass  mode  by  White,  Malmberg,  and  Driscoll^^.  For  m  >  1, 
however  the  modes  are  already  damped  by  the  collisionless  phase  mixing  just 
described.  When  external  resistors  are  added  in  the  m=2  case  as  shown  in 
Fig.  12c,  the  collisionless  damping  rate  is  decreased .  Fig.  15  shows  the 
damping  rate  versus  resistance.  Furthermore,  with  sufficient  external 
dissipation  the  otherwise  damped  mode  is  turned  in  to  a  growing  mode. 

We  have  seen,  by  examining  the  electric  field  at  the  conducting  wall, 
that  the  linear  response  to  a  short  applied  pulse  is  a  decaying  sinusoid, 
approximately  exponential.  However,  density  perturbations  at  each  radius 
within  the  plasma  persist,  indefinitely  in  the  ideal  drift  approximation. 

This  can  be  seen  by  inverting  Eq.  8 
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dn 

=  Ti-dr 


J  [w-mw.Cr)]  ®  27r 


(13) 


which,  because  of  the  simple  pole  at  w  =  mw^(r)  gives  an  undamped 
contribution  to  n^(r,t)  proportional  to  exp[-im(*;^(r ) t] .  Thus,  perturbations 
at  each  radius  continue  to  be  convected,  undamped,  at  the  local  angular 
velocity.  The  density  N^(r,w)  and  n^(r,t),  calculated  by  the  numerical 
methods  described  above,  is  shown  for  several  representative  radii  in  Fig. 

16.  It  is  seen  that  nj(r,t)  does  not  decay  (lower  panel)  but  has  different 
frequencies  (upper  panel)  at  different  radii.  However,  the  electric  field  at 
the  conductor  (and  therefore  the  charge  on  the  conductor)  decays  away  because 
it  is  caused  by  (a  weighted  average  of)  charge  density  perturbations  at 
various  radii.  This  result  is  expressed  mathematically  as  the  derivative  with 
respect  to  r  of  the  integral  in  Eq.  6.  Because  of  the  radially  dependent 
phase  factor  in  the  integrand,  contributions  from  different  radii  have 
different  frequencies  and  phase  mix  away.  The  initial  and  phase  mixed 
density  perturbations  are  illustrated  in  Figs.  IjSa  and  1^,  respectively. 


V.  FLUID  ECHOS 

The  dynamics  of  the  ideal  fluid  is  time  reversible  and  one  suspects  by 
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analogy  with  the  cyclotron  echos  and  plasma  wave  echos  that  this  system 
might  also  exhibit  echos.  The  application  of  a  second  pulse  might  partially 
unmix  the  undamped  phase-mixed  remnants  of  the  first  pulse  and  create  a 
coherent  response  at  a  later  time,  long  after  the  responses  of  both  pulses 
have  decayed  away.  The  mixing  occurs  because  of  different  angular  velocities 
so,  by  analogy  with  the  plasma  wave  echo,  we  examine  the  effects  of  two  pulse 
applied  to  the  electrode  structure  shown  in  Fig.  12.  The  connections  are 
arrranged  so  as  to  produce  potentials  with  different  angular  dependences:  m^ 
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for  the  first  pulse  and  m  for  the  second  pulse.  For  simplicity  we  take  the 

A 

time  dependence  of  the  applied  pulse  to  be  V^t^^(t)  and  V^t^^Ct-r) 
respectively.  This  form  serves  as  the  Green’s  function  (in  time)  for  other 
pulse  shapes.  We  use  the  following  notation: 

4>  =  +  <f>^(T,t)  (14a) 

n  =  n^(r)  +  n^(r,t)  +  n2(r,t)  +  ng(r,t)  (14b) 

where  the  subscript  o  refers  to  the  steady  state,  subscripts  i  and  2  refer  to 
the  first  and  second  applied  pulses,  and  subscript  3  refers  to  the  echo, 
which  results  from  the  nonlinear  interaction  between  pulses  1  and  2. 

We  treat  the  first  and  second  pulses  in  the  linear  approximation  (as  described 
above),  and  seek  the  the  nonlinear  interaction  between  the  two  pulse  to  obtain 
the  echo  response  to  second  order  in  the  amplitudes.  To  this  order  the 


continuity  equation  for  the  echo  may  be  written 


3  W  3 

dt  ^  dd 


dn. 

0 

+  3r  j — 
dr 


dr  ^  r  d6  ^  *  ~r  dT 


where  the  nonlinear  terms  arising  out  of  the  interaction  of  the  pulses  1  and  2 
have  been  put  on  the  right  hand  side.  These  terms,  which  we  denote  by 
S3(r,t),  act  as  source  terms  for  the  echo  and  are  obtained  from  the  linear 
solutions  for  pulses  1  and  2.  The  left  hand  side  is  the  usual  linearized 


continuity  equation. 


Since  we  are  using  complex  exponentials  to  describe  quantities  which 
are  real  and,  by  convention,  taking  the  real  part  (or  adding  the  complex 
conjugate  and  dividing  by  2),  some  care  is  required  with  the  quadratic  source 
term.  It  is  readily  shown  that  when  using  complex  notation  one  should  use  the 


complex  conjugate  of  the  complex  solution  for  one  of  the  pulses  and  divide  by 
2.  We  chose  to  do  this  for  pulse  1,  so  that  the  echo  has  angular  mode  number 


m  =  m  -  m  >0,  where  m  and  m  are  positive.  Since  the  source  term  s  (r,t) 
consists  of  products  of  time-dependent  quantities,  the  Laplace  transform  of 
s  is  given  by  the  convolution  of  the  transforms  of  the  individual  quantites. 

Of  the  four  terms  in  s  ,  the  first  dominates  because  n  (r,t)  contains  terms  ~ 
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exp[-ia;Q(r)t] .unlike  Vj(r,t).  Because  it  is  differentiated  with  respect  to  r, 
it  contains  terms  which  grow  with  t.  In  transform  space  the  source  term  has 
the  form 


S  (r.w)  « 

o 


n 

J 


$  (r.w-u;') 


dNjlr.w') 

dr 


do;* 

2?r 


(16) 


in  the  dominant  term  approximation. 


We  now  describe  the  method  of  solution.  Eqs.  (7)  and  (8)  must  be  solved, 

for  a  specified  density  profile,  to  obtain  $  (r.w),  N  (r,a;),  $  (r,Ci;),  and 
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N  (r.w)  for  pulses  1  and  2.  Analytic  solutions  can  only  be  obtained  for 
2 

special  profiles  so  we  resort  to  numerical  methods  and  consider  discrete 
values  of  r  and  w,  rj  =  j  Ar,  =  k  Atv  +  i(7.  In  particular  we  use  512  or 
1024  discrete  values  of  w  so  as  use  DFTs  to  obtain  the  time  functions.  For 
each  of  the  discrete  values  of  w,  we  numerically  integrate  Eqs.  (7)  and  (8) 
from  r=0  to  r=a  using  a  fourth  order  Runge-Kutta  method  so  as  to  obtain  values 
of  the  functions  at  500  radial  points  within  the  plasma.  Since  we  also  need 
dNj(r,a;)/dr  it  is  also  calculated  during  the  numerical  integration  at  the  500 

radial  points.  The  potentials  so  obtained  must  vary  as  (r)™!  and  (r)”’2  near 

the  origin  and  match  the  Laplace  transforms  of  the  applied  potentials  at  the 

cylinder,  respectively.  The  convolution,  Eq.  (16)  ,  is  then  done  numerically 

to  obtain  S  (r.,a;,).  The  equations  for  the  echo  pulse  are  similar  to  those  of 
3  j  k 

pulses  1  and  2,  except  that  instead  of  being  driven  by  potentials  applied  to 


the  conducting  cylinder,  the  source  term  is  S  : 

O 
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(17) 
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These  are  also  integrated  numerically  in  radius  for  512  discrete  frequencies 
and  inverted  with  an  DFT  to  obtain  the  time  dependence.  Fig.  18  shows  the 
response  to  pulses  1  and  2  as  well  as  the  echo  for  a  typical  profile.  Here 
m  =2,  m  =4,  and  m  =2.  In  each  case  the  time  dependent  electric  field  at  the 
wall  is  shown.  The  echo  amplitude  is  proportional  to  the  product  of  the 
amplitudes  of  pulses  1  and  2,  V  t  V  t  ,  is  generally  much  smaller  than  the 

L  1  Z  Z 

responses  to  the  two  pulses.  It  is  not  shown  to  scale  in  Fig.  18. 


The  echo  amplitude  is  expected  to  be  diminished  by  viscous  effects  as  the 
time  between  the  first  pulse  and  the  echo  lengthens.  This  is  because 
perturbations  at  nearby  radii  have  been  convected  by  substantially  different 
amounts  in  azimuth  causing  the  radial  gradients  become  very  large.  Thus  the 
echo  phenomenon  may  form  the  basis  for  studying  viscous  effects.  However,  the 
particles  primarily  responsible  for  the  echo  are  restricted  to  the  outer  radii 
of  the  plasma,  so  that  is  probably  the  only  region  which  could  be  probed. 

VI.  HIGH  FREQUENCY  MODES. 

So  far  the  discussion  has  been  limited  to  low  frequencies  where  the  drift 
approximation  is  valid.  2D  modes  near  the  cyclotron  frequency  also  exist 
whose  frequency  is  given  by  Eq.  2b  for  uniform  density  (rigid  rotor  steady 
state).  The  m=l  center  of  mass  mode  lies  below  the  cyclotron  frequency  by  the 
diocotron  frequency,  and  the  higher  modes  are  up  shifted  by  the  rotation  of 
the  plasma.  In  this  section  we  again  focus  on  nonuniform  density  and  angular 
velocity  profiles  such  as  shown  in  Fig.  7b.  Inertial  terms  can  no  longer  be 
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neglected,  and  expressions  for  the  fluid  velocity  in  a  nonuniformly  rotating 
fluid  become  much  more  complicated. 

4  5 

An  analysis  of  the  high  frequency  behavior  ,  near  the  cyclotron 
frequency,  of  a  cold  nonuniformly  rotating  pure  electron  plasma  has  been 
carried  out.  The  mode  structure  and  the  admittance  function  Y(w),  as  defined 
in  Sec.  IV,  has  been  evaluated  for  a  number  of  stable  monotonic  profiles.  The 
results  are  not  unlike  those  for  the  low  frequency  diocotron  modes  described 
in  Sec.  IV.  For  m=l  there  is  a  single  discrete  mode  lying  slightly  below  the 
cyclotron  frequency 

w  =  w  -  <u;^>/2a;  .  (19) 

c  pc 

For  each  m  greater  than  1  there  is  an  absorbt ion  band 

(m-l)w  (a)  <  (w-w  )  <  (m-l)w^(O)  ,  (20) 

0  c  0 

upshifted  from  the  cyclotron  frequency  by  the  Doppler  effect  of  rotation.  As 
with  the  low  frequency  modes,  an  alternative  point  of  view  is  that  the 
absorbtion  band  is  associated  with  continuum  of  singular  normal  modes.  A 
typical  result  for  the  cold  plasma  admittance  function  is  shown  in  Fig.  19. 
These  are  positive  energy  modes  as  demonstrated  later. 

Experimentally,  the  m=l  discrete  mode  is  found,  as  expected,  below  the 
cyclotron  frequency,  but  the  m>l  absorbtion  bands  are  not  found.  Instead,  in 
the  parameter  range  where  the  absorption  bands  are  expected  from  the  cold 
plasma  theory,  a  set  of  discrete  modes  for  each  m  greater  than  1  is  found. 

Fig.  20  shows  a  two  dimensional  scan  in  density  (horizontal  axis)  and 
magnetic  field  (vertical  axis).  The  modes  for  m>l  have  been  interpreted  as 
radially  trapped  (by  the  density  and  angular  velocity  profiles)  azimuthally 
propagating  Bernstein  modes.  Bernstein  modes  arise  out  of  finite 


temperature  (finite  Larmor  radius)  effects,  which  were  not  included  in  the 
cold  plasma  rotating  fluid  analysis.  However  an  approximate  modification  of 
the  theory  which  includes  electron  temperature  has  been  given,  ’  and  the 
temperature  required  to  match  the  experimental  observations  is  consistent  with 
other  determinations  of  plasma  temperature  (0.5  -  2.0  eV) . 

For  ra>l  the  effect  of  plasma  temperature  is  to  turn  the  continuous 
absorbtion  band  described  by  Eq.  (10)  into  a  band  in  which  Bernstein  modes 
can  propagate.  This  is  similar, in  a  general  way,  to  the  effect  of  temperature 
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on  the  modes  of  a  neutral  plasma  column:  Tonks-Dattner  modes,  which  are 
propagating  Langmuir  oscillations  trapped  by  the  radial  density  profile. 
Similarly,  in  a  magnetized  neutral  plasma  the  Buchsbaum-Hasegawa  modes^^ 
near  the  2nd  harmonic  of  the  cyclotron  frequency  arise  from  Bernstein  waves 
which  are  trapped  because  of  the  radial  density  gradient.  In  both  cases  the 
inclusion  of  temperature  raises  the  order  of  the  differential  equation  and 
results  in  additional  "thermal"  modes  in  the  region  where  there  would 
otherwise  be  a  continuum  of  singular  eigenmodes. 

In  the  remainder  of  this  section  we  concentrate  on  the  m=l  center  of 
mass  mode  whose  frequency  is  just  below  the  cyclotron  frequency.  We  describe 
a  measurement  of  spontaneous  thermal  excitation  of  this  mode  and  a  method  for 
using  this  measurment  to  determine  electron  temperature.  The  method  consists 
of  measuring  the  noise  power  delivered  to  a  low  noise  reciever  caused  by  the 
fluctuating  currents,  associated  with  thermal  excitation  of  the  m=l  mode, 
which  are  induced  on  the  sector  electrodes  surrounding  the  plasma.  The 
received  power  is  compared  with  the  power  from  a  calibrated  reference  noise 
source,  and  the  ratio  of  the  powers  is  denoted  by  p.  The  absorbtion 
coefficient.  A,  is  also  required  to  determine  the  temperature.  This 


proceedure  results  in  nondestuctive  method  of  temperature  determination  for 
pure  electron  plasmas.  In  related  measurements,  spontaneous  emission 
associated  with  the  m>l  Bernstein  modes  has  also  been  observed  and  the 
inferred  temperatures  are  in  the  same  range,  although  these  measurements  are 
subject  to  greater  shot-to-shot  irreproducibility.  Our  method  is  a  variation 
on  the  method  used  by  Stenzel^®  for  neutral  afterglow  plasmas  and  is  related 
to  the  stored  ion  calorimeter  discussed  by  Wineland  These  situations 
differ  drastically  from  the  usual  radiation  transport  regime  in  astrophysics 
and  in  fusion  devices  where  radiation  is  emitted  and  absorbed  over  many 
wavelengths.  This  radiation  also  differs  markedly  from  Larmor  radiation  from 
a  collection  of  independent  single  particles,  since  we  are  dealing  with  a 
single  collective  mode  of  the  plasma. 

A  schematic  diagram  of  the  experimental  configuration  is  shown  in  Fig. 

21.  To  increase  the  sensitivity  to  the  noise  fluctuations  the  signals  from 
all  eight  sectors  of  an  octupole  section  are  combined,  with  phasing  so  as  to 
be  sensitive  to  the  m=l  mode,  and  sent  to  a  low  noise  (NF  «  3db)  140  MHz 
reciever  whose  bandwidth  is  12  kHz.  The  cyclotron  frequency  is  set  slightly 
above  the  receiver  frequency  so  that,  as  the  plasma  density  decays,  the  m=l 
mode  frequency  sweeps  upward  through  the  receiver  frequency.  By  stepping  the 
magnetic  field,  and  thus  the  cyclotron  frequency,  by  a  small  amount  on 
successive  shots,  the  m=l  mode  resonance  can  be  made  to  sweep  through  the 
receiver  frequency  at  various  times  in  the  plasma  decay.  Emission  and 
absorbtion  measurements(using  the  r.f.  generator)  are  made  on  alternate  shots. 
This  is  illustrated  in  Fig.  22,  which  shows  both  the  noise  emission  and  the 
reflected  power  (inverted  trace)  as  a  function  of  time  for  3  representative 
different  values  of  the  magnetic  field.  At  the  beginning  and  end  of  an 
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eraission  trace,  the  receiver  is  switched  briefly  to  reference  noise  source  at 
temperature,  All  elements  of  the  system  are  50  ohm  coaxial  components 

and  the  insertion  loss  of  all  elements  must  be  accounted  for. 

Typical  experimental  results:  emission,  absorbtion,  and  radiation 
temperature,  are  shown  in  Figs.  23a-c.  These  are  obtained  by  fitting 
Lorentzian  curves  to  the  type  of  data  shown  in  Fig.  22  and  using  the  maximum 
values.  The  results  from  about  20  shots  for  the  same  parameters  are  averaged 
to  obtain  more  precise  values.  Fig.  23a  shows  the  absorption  coefficient, 

A,  the  fraction  of  the  power  absorbed  by  the  plasma  when  the  r.f.  signal 
generator  is  on.  Since  the  the  errors  in  the  absorption  data  are  small  a 
polynomial  fit  is  shown  and  it  is  used  in  determining  the  temperature.  Note 
also  that  A  is  close  to  unity!  Fig.  23b  shows  the  noise  power  received  from 
the  plasma,  relative  to  that  from  the  calibrated  noise  source,  p.  The 
radiation  temperature  of  the  plasma  shown  in  Fig.  23c  is  then  determined  from 
the  relationship  T^^^  =  ip/A)  The  temperature  is  0.5  -  0.6  eV  and 

changes  only  slightly  over  400  msec;  whereas  the  plasma  density  decays  by 
nearly  a  factor  of  10  during  this  time.  The  near  constancy  of  the  temperature 
has  been  observed  on  most  measurements,  although  the  value  of  the  temperature 
varied  between  0.5  and  3.0  eV  with  different  injection  conditions.  These 
temperatures  are  typical  of  those  obtained  by  other  methods. 

Because  the  cylinder  radius  (2.5cm)  is  small  compared  to  the  free  space 
wavelength  (215cm)  of  the  radiation,  these  results  can  be  interpreted  in 
terms  of  the  electric  circuit  model  shown  in  Fig.  24.  If  the  plasma  is  at 
temperature  T,  the  power  emitted  per  unit  bandwidth  into  the  load  (500)  is 
/cTA,  where  k  is  Boltzmann’s  constant  and  A  is  the  absorbtion  coefficient. 

The  latter  is  obtained  from  the  equivalent  circuit  as  A  =  4R^R/ |R^-i-Z(u;)  |  ^ , 


# 


where  the  impedance  function,  Z(w)  =  1/Y(w),  is  the  inverse  of  the  admittance 
function. 

One  of  the  surprising  results  is  that  the  m  =  1  mode  absorbs  energy 
since  it  is  expected,  from  elementary  considerations,  to  be  lossless. 
Furthermore,  as  illustrated  in  Fig.  23a,  the  absorption  coefficient  changes 
only  slightly  while  the  density  decays  by  nearly  an  order  of  magnitude.  It  is 
possible  to  obtain  the  frequency  shift  and  linewidth  versus  time  of  the  m=l 
mode  from  the  Lorentzian  fits  to  the  data  of  Fig.  22.  These  are  shown  in 
Fig.  25.  The  linewidth  is  10-40  kHz,  a  significant  fraction  of  the  frequency 
shift,  30-180  kHz.  Furthermore  the  linewidth  tracks  the  density  decay.  The 
frequency  width  of  the  mode  is  found  from  Aw=(dw^/dt)At,  where  At  is  the  time 
to  sweep  through  the  mode  as  determined  from  Fig.  22  and  dw^/dt  is  the  rate  of 
change  of  the  m=l  mode  frequency.  Since  the  absorption  coefficient  A  is  about 
unity,  the  plasma  is  matched  to  the  5011  system,  and  therefore  half  the 
linewidth  is  due  to  the  external  system  and  the  remaining  half  is  intrinsic  to 
the  plasma. 

An  independent  determination  of  the  line  width  of  the  m=l  cyclotron  mode 

has  been  made  by  measuring  its  decay  rate  when  excited  by  a  short  burst  at  the 

mode  frequency.  By  adding  5011  resistors  to  various  numbers  (N)  of  the 

octupole  sectors,  their  effect  on  damping  rate  can  be  found.  This  is  the  high 
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frequency  analog  of  the  m=l  and  m=2  low  frequency  resistor  experiments,  ’ 
except  that  here  the  external  dissipation  increases  the  damping,  indicating 
that  the  mode  is  a  positive  energy  mode.  The  result  is  shown  in  Fig.  26  and 
shows  clearly  that  the  damping  rate  increases  linearly  with  N.  The  minimum 
value  of  N  is  2,  from  the  source  and  the  receiver,  giving  a  decay  time  of 
about  50  nsec.  To  determine  the  intrinsic  (N=0)  damping  rate  of  the  mode  we 
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plot  the  damping  rate  versus  N  and  determine  the  intercept  corresponding  to 
N=0.  This  result  is  shown  in  Fig.  27  for  several  different  densities.  The 
most  important  conclusions  are  that  the  intrinsic  damping  rate  is  non-zero  and 
that,  for  low  densities,  the  intrinsic  damping  rate  increases  with  density. 

The  latter  is  consistent  with  A  being  approximately  independent  of  density. 
VII.  DISCUSSION. 

Experimental  results  presented  above  suggest  some  unresolved  issues 
which  may  lead  to  a  still  better  understanding  of  the  phenomena. 

For  m=l  the  cold  plasma  equations  for  a  nonneutral  plasma  predict  two 

discrete  undamped  center  of  mass  modes.  The  low  frequency  mode  is  observed  to 
2  3  5  2 

be  undamped  or  weakly  damped  ;  whereas,  the  high  frequency  mode  is  strongly 
damped.  This  origin  of  this  damping  is  as  yet  unexplained. 

For  m>l ,  the  cold  plasma  equations  predict  absorption  bands  (continuum  of 
singular  eigenmodes)  for  both  the  low  frequencies,  w  ~  and  high 
frequencies,  However,  experimentally,  radially  trapped  Bernstein 

modes  are  found  in  the  frequency  band  in  which  the  high  frequency  absorption 
band  is  supposed  to  appear.  These  require  temperature  or  finite  Larmor  radius 
(FLR)  effects.  However,  temperature  and  FLR  effects  appear  to  be  unimportant 
for  the  the  low  frequency  modes.  Why? 

For  the  collisionless  damping  of  low  frequency  disturbances  in  sheared 
flows,  such  as  those  described  in  Sec.  VI,  what  is  the  role  and  importance 
of  viscosity.  One  normally  thinks  of  viscosity  as  a  dissipative  phenomenon 
and  we  have  seen  that  external  dissipation  reduces  the  collisionless  damping 
rate.  Does  viscosity  have  the  same  effect? 

Approximate  solutions  of  the  Vlasov  equation  in  sheared  flow,  such  as 
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those  presented  by  Pearson  for  neutral  cylindrical  plasmas  and  by  Prasad  et 
al^^  for  a  nonneutral  slab  plasma,  are  needed  for  both  the  high  frequency  and 
low  frequency  regimes.  In  order  to  carry  out  such  analyses  it  is  necessary  to 
find  slowly  evolving  self-consistent  sheared  profiles  from  which  to  start. 

Plasma  temperature  appears  to  decay  much  more  slowly  than  plasma  density 
in  pure  electron  plasmas.  This  observation  needs  an  explanation.  Release  of 
energy  through  expansion  of  the  plasma  as  the  density  decays  probably  plays  a 
role  in  determining  the  temperature. 
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Advances  in  Nonneutral  Plasmas 

Brillouin  Flow.  a.  cylindrical,  b.  planar,  b.  is  shown  with  a 
slow-wave  circuit  as  used  in  Backward  Wave  Oscillators. 

Penning  traps,  a.  cylindrical  and  b.  hyperbolic.  Radial  containment 
is  by  the  axial  magnetic  field  and  axial  containment  is 
electrostatic . 

Rotation  Frequency  versus  plasma  density  in  Brillouin  Equilibrium. 
Experimental  points  are  from  Ref.  8. 

Surface  modes  of  a  cylindrical  nonneutral  plasma  with  constant 
density. 

Growth  rate  versus  thickness  for  planar  Brillouin  beam. 

Images  of  unstable  hollow  plasmas  showing  azimuthal  wavelength 
for  three  thicknesses  (0.67,  0.17,  and  0.14  cm).  Left  side: 
before  instability,  Right  side:  after  onset.  From  Ref.  33. 

Density  profiles  and  corresponding  angular  velocity  profiles, 
monotonical ly  decreasing  and  stable  (upper),  and  hollow  profiles 
which  are  unstable  (lower). 

Measured  z-averaged  density  n(r,ff)  at  t=  50,  120,  170,  and  1000 
fisec,  as  m=2  hollow  profile  instability  grows,  saturates,  and 
evolves  to  a  stable  monotonic  profile.  From  Ref.  34. 

Fluid  code  calculation  of  n(r,0,t)  at  6  successive  times  in  the 
evolution  (upper  left  to  lower  right)  of  the  m=2  hollow  profile 
instability.  Density  scale  is  shown  at  the  right. 


Fig.  10.  a.  Bromwich  contour  in  complex  w  plane  used  in  inverting  the  Laplace 
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Fig. 


Fig. 
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Fig. 


Fig. 


transform,  b.  Deformed  contour  revealing  simple  pole  in  the  lower 
half  plane. 

11.  (a)Real  and  Imaginary  Parts  of  (rE^/(^)j._jj  as  a  function  of  frequency 

normalized  to  the  central  rotation  frequency,  for  m=2.  The  vertical 

dotted  line  corresponds  to  the  m=2  resonant  frequency  w  =  .904  for 

this  profile,  (b)  Semilog  plot  of  the  envelope  of  inverse  Fourier 

transform  of  (rE  /T<i>)  ,  .  The  initial  decay  changes  from 

r  r=D 

exponential  to  algebraic.  From  Ref.  36. 

12.  (a)  Schematic  of  cylindrical  structure  for  plasma  excitation  and 
trapping,  (b)  phasing  of  first  octupole  for  exciting  an  m=2 
disturbance,  (c)  configuration  of  second  octupole  for  signal 
reception,  and  (d)  configuration  of  second  octupole  for  negative 
energy  test.  From  Ref.  36. 

13.  Decay  of  the  envelope  of  m=2  signal  at  about  400  KHz  excited  by  and 
r.f.  pulse  with  voltages  10  (lower  left),  20,  40,  60,  80,  120,  160, 
200,  250,  300,  400,  500,  600,  700,  and  1000  (uppermost)  mV,  plotted 
on  a  semi-log  scale  with  vertical  scale.  From  Ref.  36. 

14.  2D  Fluid  Code  Calculation  of  Damping  and  Trapping.  Plotted  is  the 
charge  induced  on  one  octant  of  the  cylinder,  and  in  each 
successive  trace  the  excitation  amplitude  is  doubled. 

15.  Decrease  in  Decay  Rate  of  the  m=2  response  versus  resistance  when  a 
a  resistance  R  is  connected  to  the  wall  electrode.  The  line  is  a 
least  squares  fit  to  the  experimental  data.  From  Ref.  36. 

16.  Density  perturbations  at  (r/a)  =  0.93,  0.87,  0.81,  0.75,  0.69,  and 
0.63  (traces  a-e).  Upper  figure  is  magnitude  of  Laplace  transform 
of  the  density  along  the  contour  in  Fig.  10a  and  the  lower  figure  is 


a  semilog  plot  of  density  versus  time. 

Fig.  17.  Internal  in=2  perturbation  of  cylindrical  nonneutral  plasma,  (a) 

Density  countours  when  perturbations  at  all  radii  are  in  phase  and 
(b)  after  phase  mixing. 

Fig.  18.  Time  dependence  of  the  response  to  two  applied  pulses,  m=2  at  t=0  and 
m=4  at  t=r,  followed  an  m=2  echo  at  t=2r.  Upper  curve:  Waveforms  of 
the  electric  field  at  the  wall,  lower  curve:  semilog  plot  of  the 
magnitudes  of  waveforms.  The  echo  amplitude  has  been  increased  for 
clarity. 

Fig.  19.  Real  and  Imaginary  parts  of  the  admittance  function  versus  frequency 
for  high  frequency  modes  (T=0).  When  Re{Y(w)>0  absorption  occurs 
due  to  an  upper  hybrid  resonance  somewhere  within  the  plasma. 

Due  to  rotation,  absorption  bands  are  upshifted  by  rotation.  From 
Ref.  41. 

Fig.  20.  Experimental  Scan  of  k  =  0  modes  near  the  cyclotron  frequency, 
showing  a  single  m=l  center  of  mass  mode  below  the  cyclotron 
frequency,  and  bands  of  m=2  and  3  Bernstein  modes  above. 

Fig.  21.  Schematic  for  measurement  of  radiation  temperature  from  the  m=l 

mode.  Emission  and  Absorption  are  measured  on  alternate  shots  and 
before  and  after  the  emission  measurement  a  comparison  with  a 
calibrated  noise  source  is  made. 

Fig.  22.  Einission  and  absorbtion  traces  versus  time.  As  the  density  decays 
the  resonance  sweeps  through  the  receiver  frequency  and  by  changing 
the  magnetic  field  slightly  the  resonance  can  be  positioned  at 
various  times.  Lowest  trace  is  semilog  plot  of  density  with  each  tic 
mark  indicating  a  factor  of  10. 
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Fig. 

Fig. 


Fig. 


Fig. 


Fig. 


23.  Emission,  Absorbtion,  Temperature  versus  time  obtained  from  35 
traces  similar  to  those  shown  in  Fig.  22. 

24.  Model  for  emission  and  aborbtion.  Left,  electrical  schematic  of  the 
of  the  experiment.  Right,  equivalent  circuit  with  the  plasma 
replaced  by  equivalent  noise  source  and  its  impedance  function. 

25.  Experimentally  determined  frequency  downshift  and  linewidth  versus 
time  form  the  m=l  mode,  as  determined  from  traces  similar  to  those  of 
of  Fig.  22. 

26.  Semilog  plot  of  the  free  decay  of  the  m=l  mode  amplitude  versus  time 
when  excited  by  a  short  resonant  burst.  Traces  are  for  2  to  16 
sectors  terminated  in  5011,  indicating  that  additional  traces 
increase  the  decay  rate. 

27.  Measured  free  decay  rate  of  the  m=l  mode  versus  the  number  of 
sectors  with  resistors.  The  non-zero  intercept  indicates  an 
intrinsic  damping  of  the  m=l  mode. 
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THEORY  OF  CYCLOTRON  RESONANCE  IN  A  CYLINDRICAL  NON-NEUTRAL  PLASMA 


Roy  W.  Gould 
Applied  Physics 

California  Institute  of  Technology ,  Pasadena ,  CA  91125 


Resonances  near  the  cyclotron  frequency  of  a  cylindrical  non-neutral  plasma 
column  with  radial  density  and  angular  velocity  profiles  are  studied.  Cold 
plasma  and  warm  plasma  models  are  used  to  consider  the  various  multipole  modes 
for  =  0.  For  m=l  only  a  single  resonance,  a  center  of  mass  mode,  is 
found.  Its  frequency  is  downshifted  from  the  cyclotron  frequency  by  an  amount 
equal  to  the  frequency  of  the  low  frequency  diocotron  mode.  For  each  m 
greater  than  1  the  cold  plasma  model  gives  a  continuous  absorption  band  with 
spatially  localized  absorption  at  each  frequency,  corresponding  to  a  continuum 
of  singular  modes.  In  the  warm  plasma  model,  the  continuuo  go  over  to  sets  of 
discrete  radially-trapped  azimuthal ly-propagating  Bernstein  modes.  For  both 
cases  the  plasma  rotation  Doppler-shifts  the  azimuthally  propagating  modes  up 
in  frequency.  The  angular  velocity  profile,  together  with  plasma  temperature, 
determine  the  spacing  of  the  Bernstein  modes  within  a  band. 
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I.  INTRODUCTION 


In  this  paper  the  inodes  of  a  cylindrical  nonneutral  plasma  column  in  the 
vicinity  of  the  cyclotron  frequency  a;^=eB/m^  are  studied  theoretically.  This 

work  was  done  in  support  of  recent  experiments.^  Plasma  rotation  plays  a  vital 

role  in  determining  the  mode  structure  near  the  cyclotron  frequency  of 

nonneutral  plasmas  and  is  responsible  for  marked  differences  from  otherwise 

similar  phenomena  in  neutral  plasmas.  Nonneutral  plasmas  have  been  studied 

2  3 

extensively  in  recent  years  and  there  are  two  good  summaries  of  this  work. 

However  there  has  been  little  attention  to  high  frequency  phenomena.  A 

3  4 

notable  exception  is  the  early  theoretical  work  of  Krall  and  Davidson  ’  in 
which  they  show  how  to  carry  over  certain  Vlasov  results  from  cylindrical 
neutral  plasmas  to  cylindrical  nonneutral  plasmas.  However,  their  results  are 
limited  to  the  rigid  rotor  equilibrium  state,  a  situation  not  always  obtained 
in  laboratory  plasmas  because  of  the  time  required  to  reach  this  state 
can  exceed  the  plasma  lifetime.  Recently  Dubin  has  given  a  theory  of  the 
modes  of  cold  rigid  rotor  spheroidal  nonneutral  plasma,  which  includes  high 
frequency  modes.  Pearson^  gave  a  finite  Larmor  radius  treatment  of  the 
Bernstein^  waves  of  a  slab  and  cylindrical  neutral  plasmas  with  a  spatially 
varying  density  profile  but  rotation  was  not  included  in  the  cylindrical  case. 
More  recently  Prasad,  Morales,  and  Fried  have  have  given  a  small  Larmor 
radius  theory  of  modes  of  a  nonneutral  slab  plasma  which  has  sheared  flow. 
However,  because  of  the  slab  geometry,  their  results  cannot  be  compared 
directly  with  experiments. 

In  this  paper  we  are  particularly  interested  in  the  modes  of  a 
non— equi 1 ibr ium  plasma  in  which  both  plasma  density  and  angular  velocity  are 
functions  of  radius.  While  they  are  expected  to  evolve  toward  thermal 


equilibrium,  this  evolution  takes  place  on  a  sufficiently  long  time  scale  that 
the  plasma  may  be  considered  to  be  in  a  steady  state  for  the  purposes  of 
describing  the  phenomena  discussed  here.  For  simplicity  we  examine  only  the 
k^=0  modes  of  a  long  column,  but  we  consider  the  various  cylindrical 
harmonics,  m?K).  First  we  discuss  the  cold  plasma  case  which  leads  to  a 
continuous  spectrum  of  modes  for  each  m.  We  then  give  an  approximate 
treatment  of  the  warm  plasma  modes  (first  order  in  temperature),  which  leads 
to  azimuthally  propagating  and  radially  trapped  Bernstein  modes,  reminiscent 

9 

of  the  Buchsbaura-Hasegawa  modes  which  occur  near  in  a  neutral  plasma 

which  has  a  radial  density  profile  but  no  rotation^°. 

It  is  instructive  as  a  preliminary  to  recall  the  behavior  of  a  cold 

cylindrical  neutral  plasma  whose  density  depends  upon  radius  but  which  is  not 

rotating.  When  an  oscillating  electric  field  is  applied  externally,  the 

plasma  exhibits  a  spatially  localized  resonant  response  at  the  radius  where 
2  2  2 

Ki(r)=l-i»;  (r)/(u;  -w  ) ,  the  perpendicular  dielectric  function,  vanishes, 
p  c 

Because  of  the  radial  dependence  of  density,  there  is  a  range  of  frequencies 
for  which  this  condition  can  be  satisfied,  corresponding  to  the  range  of  upper 
hybrid  frequencies  within  the  plasma.  This  leads  to  a  continuum  of  singular 
eigenmodes,  similar  to  those  found  in  other  problems  ’  ,  and  an  upper  hybrid 

continuum  absorption  band. 

A  non-neutral  plasma  column  is  necessarily  in  rotation  in  order  to  be  in 
a  steady  state.  This  rotation  gives  rise  to  a  Doppler  shift  of  an  applied 
frequency,  and  the  upper  hybrid  resonance,  as  seen  in  the  laboratory  frame  is 
upshifted.  As  we  shall  see  later,  the  centrifugal  and  Coriolis  forces 
associated  with  the  rotation  also  give  rise  to  a  downshift  of  the  single 
particle  gyration  frequency.  For  rigid  rotor  equilibria  the  downshifted 


frequency  is  called  the  vortex  frequency^.  For  non-rigid  rotor  steady  states, 
the  single  particle  downshift  and  the  Doppler  upshift  both  depend  upon  radius 
and  this  complicates  the  situation  considerably. 

In  the  remainder  of  the  paper  we  study  the  response  of  the  plasma  to 
applied  multipolar  fields  produced  by  the  geometry  of  Figure  1.  In  Sections 
II-VI  we  discuss  the  cold  plasma  behavior  of  electrons  and  the  absorbtion  band 
which  arises  out  of  the  continuum  of  singular  normal  modes.  In  Section  VII  we 
show,  approximately,  how  the  continuum  goes  over  to  a  set  of  radially-trapped 
azimuthal ly-propagating  Bernstein  modes  in  a  warm  plasma,  where  finite  Larmor 
radius  effects  are  important.  In  Section  VIII  we  extend  the  results  to 
multiple  ion  species.  Throughout  the  paper  it  is  assumed  that  is  a 

small  parameter,  and  results  are  obtained  only  to  first  order  in  this 
parameter.  We  also  assume  that  the  dimensions  are  small  compared  with  a  free 
space  wavelength  so  that  retardation  effects  can  be  neglected. 

II.  COLD  PLASMA  EQUATIONS. 

We  consider  a  cylindrical  non-neutral  plasma  which  is  contained  by  a 
static  axial  magnetic  field  Bq  and  which  has  a  radial  density  profile  n^Cr) 
which  is  assumed  to  be  known.  The  density  profile,  ngCr),  need  not  be  an 
equilibrium  profile  in  the  thermodynamic  sense,  since  a  variety  of 
non-equilibrium  profiles  can  be  produced  which  evolve  on  a  time  scale  very 
long  compared  to  the  time  scale  of  cyclotron  resonance  observation  times. 
Associated  with  a  given  density  profile,  n^Cr),  is  a  radial  electric  field 
profile,  E^Cr),  and  an  angular  velocity  profile,  w^Cr),  which  we  describe  in 
Section  III. 

The  cold  fluid  equations  of  motion  and  continuity  in  Eulerian  form  for 


cylindrical  geometry  and  no  z  variation  are: 


^  Idi  +  u  du  Idr  +  (u^/rldu  IdO  -  u^/r  =  -(e/m  )[E  +  Uo  B_],  [1] 

r  rr  u  x  u  ert/u 

du^/dt  +  \i_dagldr  +  (u^/r)du^/50  +  u^u^/r  =  -(e/m^)[E^-  B^],  [2] 

dvildt  +  (l/r)d(rnu^)/dr  +  (l/r)5(nu^)/9fl  =  0.  [3] 

Maxwell  equations  for  the  quasi-static  electric  field  are: 

(l/r)d(rE^)/dr  -  {llr)dE^lde  =  0,  [4] 

(l/r)d(rE^)/ar  +  (l/r)dE^/d^  =  -  ne/e^.  [5] 

III.  THE  STEADY  STATE. 


In  the  steady  state,  E^  and  u^  vanish  and  the  remaining  quantities,  E^, 
and  n  are  independent  of  6  and  t.  The  steady  state  quantities  obey  the 
equations : 


d(rEQ^)/dr  =  -  rn^e/e^  , 


[6] 

[7] 


where  a;Q(r)=UQ^(r)/r  is  the  angular  velocity  at  radius  r.  The  subscripts  0 

refer  to  the  steady  state  values.  Since  we  are  concerned  with  the  low  density 
2  2 

limit,  (ci;p/ci;^)«l,  the  rotational  angular  velocity,  is  everywhere  small 

2 

compared  with  the  cyclotron  frequency  and  the  centrifugal  force  term  rw^  of 
[6]  may  be  neglected.  This  corresponds  to  being  well  below  the  Brillouin 
density  limit.  Eqs.  6  and  7  then  lead  to  the  usual  low  density  relationship 
between  angular  velocity  and  density: 


(l/r)d[r^a;„(r)]/dr  =  w^(r)/a;  .  [8] 

0  pc 

We  find  it  convenient  to  normalize  the  radial-dependent  plasma  frequency  and 
rotational  angular  velocity  to  their  central  values  and  introduce  the 
functions  f(r)  and  g(r)  to  express  the  normalized  density  and  angular  velocity 


functions.  We  let 
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J^(t)  =  w^(0)  f(r) 
P  P 

WoCr)  =  Wq(0)  g(r) 


[9] 

[10] 


where,  by  definition,  f(0)  =  g(0)  =  1.  It  follows  from  Eq.  8  that 

the  central  rotation  frequency  is  ^^(O)  =  t.;p(0)/2w^.  With  these  definitions 

Eq.  8  takes  the  form 

d[r\(r)]/dr  =  2  r  f(r).  [8a] 

A  simple  illustration  of  a  monotonically  decreasing  density  profile  is  the 
parabolic  density  profile  f (r)=l-(r/a)^ ,  for  r<a.  This  profile,  together  with 
the  corresponding  angular  velocity  profile,  g(r)=l-.5(r/a)  .  will  be  used  in 
later  discussions,  although  the  general  approach  is  not  limited  to  this 
example.  g(r)  is  obtained  by  integrating  Eq.  8a  and  both  f(r)  and  g(r)  are 
shown  in  Fig.  2. 

IV.  SMALL  PERTURBATIONS. 

The  linearized  forms  of  Eqs.  1-3  for  small  perturbations  (subscript  1) 
when  one  assumes  exp(im0-iwt)  dependence  are: 


-i[w-ma;Q]Uj^^+  [u;^-2u;q]Uj^^  =  -(e/m^)Ej^^ 
-[w^-2a;Q-rda;Q/dr]Uj^  -i[c.;-ma;jj]u^^  =  -(e/m^lEj^^ 
-iCw-mw^lnj  +  (l/r)d(rnQUjj.)/dr  +  (im/rln^u^^  =  0 
which  are  readily  solved  for  u^^  and  Uj^^  in  terms  of  Ej^^  and 

-i[‘^“"“^o]Eir  ~ 


[11] 

[12] 

[13] 


u 


Ir 


=  -  (r  ) 


”e  -[w-mwp]'^  +  [u;^-2wQ][a;^-2a;Q-rdc.;Q/dr] 


u..  =  -  (r  ) 


'•le 


[w^-2a;Q-rdwQ/dr]Ej^  -  i[c<;-mwQ]Ej^^ 


[14] 


[15] 


e  -[w-mw^l  +  [a;^-2a;Q][a;^-2wQ-rda;Q/dr] 


Now  a;'=w-mWQ(r)  is  the  local  Doppler  shifted  frequency  seen  by  a  fluid  element 
drifting  around  the  axis  at  radius  r.  The  other  term  in  the  denominator  of 
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Eqs.  14  and  15  can  be  identified  as  the  square  of  the  effective  single 
particle  gyration  frequency 

w'  =  [  (w^-2Wq)  (w^-2wQ-rdc.;Q/dr)  [16] 


This  frequency  is  downshifted  slightly  from  w^=eBQ/m^  by  the  centrifugal  and 
Coriolis  forces  and,  in  general,  the  downshift  is  also  dependent  upon  radius. 
However,  in  the  special  case  when  the  plasma  density  is  independent  of  radius, 
Wq  is  also  independent  of  radius  and  one  recovers  the  rigid  rotor  vortex 
frequency,  a;'=a;  which  is  also  independent  of  radius. 

C  C  V/ 

Since  we  are  primarily  interested  in  frequencies  close  to  the  cyclotron 
frequency  we  can  simplify  Eqs.  14  and  15.  It  is  useful  to  introduce  the 
normalized  frequency  A  through  the  definition 

w  =  ^c  +  (Wp(0)/2a;^)  A,  [17] 

where  A  is  the  frequency  difference  normalized  to  the  central  rotation 

2  2 

frequency.  Then  to  first  order  in  w  (0)/u;  ,  Eqs.  14  and  15  can  be  written 

p  c 


"0®“lr  " 


-iu;c^of(^-)[Elr-^Eig] 
[ A-mg+2g+ . 5rdg/dr ] 


[18] 


or,  more  simply,  as  =  0  and 

-2i(jj  f(r)  E_ 

n  eu_  =  - - -  [19] 

[A-mg+g+f ] 

We  have  defined  Ej.  =  E^^^  i  iE^^^  and  Uj_  =  Uj^^  ±  iu^^^,  and  used  dg/dr  from 
Eq.  [8a].  E^  and  u^  denote  amplitudes  of  the  two  circularly  polarized 
components  of  the  electric  field  and  the  velocity.  From  this  result  we  see 
that  close  to  cyclotron  resonance  only  one  circularly  polarized  component  is 
effective  in  exciting  the  motion,  namely  that  which  matches  the  gyration  sense 
of  the  particles.  The  velocity  which  it  produces  is  circularly  polarized  with 


the  same  sense.  Linear  combinations  of  the  two  Maxwell  equations,  Eqs.  4  and 
5  give 

(l/r)d(rE^) /dr  i  (m/r)E^  =  -njC/t^  [20] 

and  the  continuity  equation,  Eq.  13,  together  with  Eq.  18  gives 

,  ,  rfE  fE 

,  _  1.  u _ -  m  _ -  ro■l^ 

”l®'^0  ”  r  dr  [A-mg+g+f]  r  [A-mg+g+f]  ’ 

where  w-mw^  w  has  been  used.  Defining  a  dielectric  function  for  right 

circularly  polarized  fields  as  K(r)  =  1  -  f /[A-mg+g+f ] ,  the  quantity  D_  =  K  E 

obeys  a  particularly  simple  equation: 

-  ^(rD  )  -  -  D  =  0  .  [22] 

r  dr  -  r  - 

D_  is  the  displacement  associated  with  the  resonantly  polarized  component. 

V.  OOLD  PLASMA  SOLUTIONS. 

111“  1 

The  solution  of  Eq.  22  is  simply  D_(r)  =  A  r  ,  where  A  is  an  arbitrary 
constant,  and  the  corresponding  electric  field  is 

E_(r)  =  A  r“"^/  K(r).  [23] 

For  m^O  A  must  be  zero  in  order  for  the  field  to  be  zero  at  r=0.  Thus  E_  is 
zero  for  m^O.  Only  applied  fields  with  positive  m  have  the  correct  circular 
polarization  to  elicit  a  response  from  a  pure  electron  plasma.  Thus  we  need 
only  consider  positive  m.  For  m=+l,  D_  =  A  is  independent  of  radius.  This 
corresponds  to  a  center  of  mass  mode. 


Having  obtained  the  solution  for  E_  we  substitute  it  into  [21]  to  obtain 
the  perturbed  density  into  Eq.  20  to  get  the  equation  for  E^: 


(rE^)  .  S  1])  .  i  ij .  [24] 


The  right  side  of  Eq  24.  gives  the  perturbed  density  in  the  plasma  and  we  can 
use  this  expression  to  obtain  the  radial  dependence  of  the  perturbed  density 
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for  various  applied  frequencies.  The  perturbed  density  exhibits  a  singularity 
at  the  radius  where  K(r)  vanishes.  This  condition  corresponds  to  a  localized 
upper  hybrid  resonance  condition,  when  the  various  frequency  shifts 
incorporated  into  K(r)  are  taken  into  account.  To  solve  for  we  note  that 

the  substitution  =  ^Kr)/r^™^^^  leads  to  an  integrable  equation  for  whose 
solution  is 


ft' [k(^-  *] 


[25] 


The  ratio  of  to  E_  at  a  given  radius  is  then  given  by 

E.(r) 


R  (r)  =  ''  =  Kill  r  1  J  ,2m,  , 

m^  E_(r)  2m  J  dr'[K(r')  ^ 


[26] 


At  a  radius  beyond  which  the  plasma  density  vanishes,  the  integral  is 

just  a  constant,  K  =  1  and  E^/E_  falls  of  as  r"^™.  Eq.  26  can  be 
integrated  by  parts  to  give 


m  2m 

r 


J„  [k(FT-T  ■ 


[27] 


for  r  >  definition  of  K(rM  this  can  also  be  written 


^ 

m  2m 

r 


Jn  U-(m-l)g(r')J’'  ’ 


[27a] 


Kjjj  is  a  measure  of  the  mth  multipole  moment  induced  in  the  plasma  column  by 

the  applied  field.  For  the  illustrative  profile  of  Fig  2  r  =  a 

’  max 

VI.  ADMITTANCE  FUNCTION. 


We  now  relate  this  result,  specifically  the  calculation  of  R  (r)  to  the 

m 

measurement  geometry  depicted  in  Fig.  1.  If  a  voltage  V  is  applied  to 

the  segmented  cylinder  shown,  the  current  I  flowing  to  the  electrode  (as 

displacement  current)  will  be 
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I  =  -iwc^  /  b  A6 


-0. 


[28] 


It  is  then  straightforward  to  show  that  the  observed  admittance  Y  =  I/V  is 


related  to  the  R  (b)  defined  above  as  follows: 
m 


_  R„(b)+1 

y(A)  =  -  i.^.^b  I 


[29] 


m 


m 


where  the  C  =  (4m/7r)[sin(m^-)/m]  are  geometrical  weighting  factors  for  the 
m  u 

various  multipole  contributions  which  depend  only  upon  m  and  the  angle  6^. 

falls  off  like  1/m  and  for  ^^=7r/2  only  the  odd  C  are  non-zero.  The  fraction 

0  m 

(R  +1)/(R  -1)  is  just  proportional  to  the  mth  multipole  component  of  E 
at  the  electrode,  and  from  Eq.  29  it  is  evident  that  we  should  examine  the 
behavior  of  that  factor  for  various  m.  Because  of  the  factor  (r'/b)^™  in  Eq. 
26  the  higher  multipole  contributions  fall  off  as  (a/b)^™,  where  we  have  set 
*^max  =  a.  Thus  the  dipole  contribution  (m=l)  will  dominate.  If  we  subtract 
the  vacuum  contribution  to  the  admittance,  i.e.  that  part  which  exists  in  the 
absence  of  the  plasma,  the  appropriate  factor  in  Eq.  29  becomes  2R^/(R^-1). 

A.  Plasma  Density  Independent  of  Radius.  For  plasma  of  uniform  density 
up  to  radius  a,  f  =  g  =  1  for  r  <  a  and  Eq.  27a  is  particularly  easy  to 
evaluate  and  gives 


yb)  =  -  (a/b)2"‘ 


SO  that 


R  (b)+l 
m _ 

R„(b)-1 

m 


A-m+l-(a/b) 


A-m+l+(a/b) 


2m 


2m 


[30] 


[31] 


Eq.  31  has  a  simple  pole  at 


A  =  m-l-Ca/b)^®  =  m-2+l-(a/b)^® 


[32] 


with  a  residue  -2(a/b)^™.  Thus  there  is  a  single  discrete  mode  for  each  m, 

with  the  coupling  decreasing  with  m,  in  accord  with  earlier  work  .  Each  term 
in  the  second  form  of  Eq.  32  can  be  identified  with  a  frequency  shift  of  a 
particular  origin.  The  first  is  a  Doppler  shift  due  to  plasma  rotation.  The 
second  is  a  downshift  of  the  single  particle  gyration  frequency.  The  third  is 
a  plasma  upshift  giving  rise  to  the  upper  hybrid  frequency.  The  last  is  a 
downshift  caused  by  the  image  charges  in  the  conductng  walls.  For  (b/a)"^  >>1 
the  influence  of  the  walls  is  negligible  and  this  shift  may  be  neglected.  In 
a  neutral  plasma,  the  first  two  shifts  are  absent.  Note  that  the  m=l 
resonance  lies  below  the  free  particle  cyclotron  frequency  and  m=2  and  higher 
resonances  lie  above. 

B.  Plasma  Density  Dependent  Upon  Radius,  m=l .  When  the  equilibrium 
density,  and  therefore  the  rotational  angular  velocity,  described  by  f(r)  and 
g(r)  respectively,  depends  on  radius  the  situation  is  more  complicated,  except 
for  m=l .  For  m=l,  the  radial  dependence  in  the  denominator  of  Eq.  27a 
disappears  and  Rj^(b)  is  readily  found  to  be 

and  Rj+1  A-<f>  [33a] 

Rj^-1  A+<f> 

where  <f>  is  the  average  of  f  over  the  entire  cross-section.  Thus  the  m=l 
resonance  occurs  at  A  =  -<f>,  i.e  it  is  downshifted  and  the  downshift  depends 
soley  on  the  density  averaged  over  the  cross  section.  This  downshift  is 

exactly  equal  to  the  frequency  of  the  low  frequency  m=l  diocotron  mode  .  In 
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Both  the  m=l  cyclotron  mode  and  the  m=l  diocotron  mode  are  center  of  mass 

modes.  The  downshift  of  the  cyclotron  mode  and  the  frequency  of  the 

diocotron  mode  both  arise  from  image  charges  in  the  conducting  wall,  as 

evidenced  by  the  fact  that  as  a/b  tends  to  zero  (wall  far  away)  both  vanish. 

Wall  charges  produce  a  downshift  because  the  force  which  they  cause  on  the 

charged  particles  is  antiparallel  to  the  X  *  ®  force,  thus  reducing  slightly 

the  frequency  of  circular  motion.  The  sum  of  the  frequencies  of  the 

downshifted  cyclotron  mode  and  low  frequency  diocotron  mode  is  equal  to 

w  =eB  /m  . 
c  0  e 

C.  Plasma  Density  Dependent  Upon  Radius,  m>l.  The  behaviour  of  the  m>l 
modes  is  not  as  simple  as  for  m=l  when  the  density  depends  upon  radius.  Then 
the  denominator  of  Eq.  27a  will  vanish  at  at  some  radius  where  A=(m-l)g(r). 
This  is  just  the  condition  at  which  the  dielectric  function  K  vanishes  and 
therefore  E  becomes  very  large.  This  condition  is  satisfied  only  at  a 
particular  radius,  which  depends  upon  the  frequency  A,  so  there  is  a  localized 
resonant  layer,  the  upper-hybrid  layer  for  a  non-neutral  plasma.  Denoting  the 
minimum  and  maximum  normalized  rotation  frequencies  where  f?^0,  then  an  upper 
hybrid  resonant  layer  occurs  in  the  frequency  band 


(m-l)g  .  <  A  <(m-l)g 

^min  ^max 


[34] 


One  can  look  at  this  another  way.  In  the  frequency  band  just  described, 

there  exists  a  continuum  of  singular  normal  modes,  with  singularity  occuring 

at  the  radius  r  where  (m-l)g(r  )=A.  Based  on  the  previous  analysis,  a 
s  s 

prescription  for  finding  the  field  of  a  singular  normal  modes  is  as  follows. 
For  a  particular  m,  pick  a  frequency  A  in  the  proper  frequency  band.  Find 
E  (r)  from  Eq.  23  by  interpreting  1/K(r)  =  P{l/K(r)}+/i5(r-rg)  where  P{  } 


stands  for  the  principal  value  and  p  is  a  constant  to  be  determined. 


Substitute  this  result  into  Eq.  25  to  find  E^(r).  Finally,  determine  the 
constant  fi  by  requiring  that  2iE^(b)=E^(b)-E_(b)=0,  since  r=b  is  a  conducting 
boundary. 

For  the  evaluation  of  the  which  appear  in  the  admittance  function  Y, 

we  can  regard  Eq.  29  as  resulting  from  the  analysis  of  an  initial  value 
problem  in  which  A  has  a  positive  imaginary  part,  thus  resolving  the  question 
of  how  to  handle  the  singularity  at  r=r  .  Alternatively,  we  could  have 

S 

assumed  the  existence  of  a  small  collision  frequency  in  the  earlier  analysis 
and  this  would  have  the  effect  of  replacing  A  by  X+iu.  For  the  purpose  of 
evaluating  and  displaying  the  various  multipole  contributions  to  the 
admittance  function,  Eq.  29,  we  follow  the  latter  approach,  taking  the  limit 
of  very  small  i/.  For  the  parabolic  profiles  previously  discussed,  the 
expression  for  Rjjj(b)  can  be  integrated  to  give 

m-2 

E„(b)=-2.(§)2"[  I  -  A-‘  -  A"-'(A-l).n(^)  ]  [35] 

k=0 


where  A  =  -2(A-m+l)/(m-l) .  When  0<A<1,  the  logarithm  is  complex  and  this 

leads  to  a  complex  admittance  function  representing  absorbtion  associated  with 

a  resonant  layer  in  the  plasma.  The  frequency  band  in  which  this  occurs  is 

just  that  given  by  Ekj.  34  with  g  .  =.5  and  g  =1.  In  particular  for  m=2, 

min  max 


R2(b)  =  -2 (|) ^ [  ( 4A-3 ) -4 (1  -A )  ( 2 A-1 ) ^n ( - ^ ^ .  [35a] 

In  Fig.  3  we  display  the  real  and  imaginary  parts  of  (R2(b)+l)/(R2(b)-l) , 
the  main  factor  in  the  admittance  Y  as  a  function  of  A.  We  also  display  the 
magnitude  of  the  perturbed  density  n^^Cr),  obtained  from  Eq.  24,  for  selected 
A.  This  clearly  displays  both  the  absorbtion  band  in  which  Y  has  a  real  part. 


and  the  existence  of  a  localized  upper  hybrid  resonance  layer  and  the 
dependence  of  its  position  on  A. 

VII.  RESONANCES  OF  THE  WARM  PLASMA. 

The  previous  discussion  assumed  zero  temperature  electrons  and  this  led 
to  singular  electric  fields  at  the  radius  where  K(r)  vanishes.  Under  these 
conditions,  finite  Larmor  radius  effects  are  important  and  lead  to  significant 
modifications  of  the  behavior  described  above.  Singular  eigenmodes  are  no 
longer  possible  and  Bernstein  modes,  propagating  across  the  magnetic  field, 
become  possible.  We  now  give  an  approximate  treatment  of  consequences  of 
finite  Larmor  radius  effects  by  modifying  the  dielectric  function  K(r) 
obtained  earlier  to  include  finite  Larmor  radius  effects,  to  first  order  in 
the  plasma  temperature.  When  the  disturbance  takes  the  form  of  a  plane  wave 
with  wave  number  k,  the  effect  of  finite  Larmor  radius  of  the  particles  is  to 
reduce  the  effectiveness  of  the  electric  field  in  producing  plasma  currents  by 
9  9  — 2  2 

a  factor  4  l^Ck  p  /2)  e  /(kp)  where  p  is  the  mean  square  Larmor 

2 

radius,  averaged  over  a  Maxwellian  distribution.  To  first  order  in  (kp)  ,  i.e 
to  first  order  in  the  temperature,  this  factor  is  approximately  1  -  .5(kp)  . 
This  approximation  is  useful  so  long  as  the  Larmor  radius  is  small  compared 
with  the  wavelength.  Since  the  particles  average  the  electric  field  over  the 
orbit,  the  dielectric  function  is  no  longer  local.  Neighboring  regions  of  the 
plasma  are  coupled  together  and  the  continuum  of  singular  normal  modes  is 
replaced  by  a  discrete  set  of  Bernstein  modes.  We  apply  this  result  to  a 
situation  where  the  waves  are  not  strictly  plane  waves  because  the  plasma 
properties  vary  slowly  with  position.  This  approximation  is  useful  so  long  as 
the  wavelength  of  the  wave  is  small  compared  with  the  distance  over  which  the 


plasma  properties  vary.  This  leads  to  the  following  modification  of  the 
dielectric  function 


K(r)  ^  1  -  ' 


[36] 


where  f  and  g  have  their  previous  meanings  and  are  functions  of  radius. 

However,  it  should  be  pointed  out  that  while  we  have  included  temperature  to 

first  order  in  the  form  of  finite  FLR  effects  in  the  wave  dynamics,  we  have 

not  included  the  effect  of  temperature  in  the  form  of  the  diamagnetic 

contribution  to  the  angular  velocity.  This  will  modify  slightly  the 

relationship  between  f  and  g.  As  before,  Eq.  36.  is  valid  only  around  the 

cyclotron  frequency.  We  now  apply  this  result  to  the  non-plane  wave  fields 

replacing  k  by  -  V  ,  the  Laplace  operator.  This  assumes  that  the  density  and 

angular  velocity  scale  lengths  are  large  compared  to  the  wavelength  of  the 

waves.  This  leads  to  an  approximate  wave  equation  for  the  Bernstein  modes 

which  will  determine  the  eigenmodes  of  the  system, 

2 


where 


r 

_  2[  (m-l)g(r)-A] 
p^f(r) 


[37] 

[38] 


k(r)  could  be  interpreted  as  the  local  wave  number  and  it  obviously  depends 

2 

upon  radius  through  g(r)  and  f(r).  When  k  (r)  >  0  the  waves  are  propagating 
and  when  k  (r)  <  0  the  waves  are  evanescent,  with  k=0  defining  the  turning 
point . 

To  proceed  further,  we  need  to  make  an  assumption  about  the  density 
profile  function  f(r).  According  to  Eq.  8a  this  also  determines  the  angular 
velocity  profile  function  g(r).  For  simplicity  we  use  the  parabolic  density 
and  angular  velocity  profiles  illustrated  in  Fig.  2.  We  note  that  this  choice 


of  f(r)  is  very  close  to  a  diffusion  profile,  f(r)  =  J^(2.405  r/a)  if  a  -  b, 
the  radius  of  the  wall.  For  (m-l)g(r)  >  A  the  Bernstein  waves  are  propagating 
and  for  (m-l)g(r)  <  A  the  waves  are  evanescent.  For  m  >  1  and  A  a  little  less 
than  (m-1)  the  waves  are  trapped  in  the  inner  part  of  the  plasma  column  and 

9 

evanscent  in  the  outer  part  of  the  column  (like  the  Buchsbaum-Hasegawa  modes 

1 3 

near  2a;  and  opposite  from  the  Tonks-Dattner  modes  near  a;  ,  which  are 
c  P 

trapped  in  the  outer  region  of  the  plasma.  The  lowest  modes  are  trapped  quite 
close  to  the  center  where  f(r)  w  1.  If  we  make  this  approximation  in  Eq.  38 
then  Eq.  37  becomes  the  two-dimensional  harmonic  oscillator  equation  for  which 
the  eigenvalues  are  readily  shown  to  be 

A,  =  (m-1)  -  (2^+m-l)  ^  ,  [40] 

where  ^  is  the  radial  mode  number. 

Alternatively,  Eq.  37  can  be  integrated  numerically  to  obtain  the 

m—l 

eigenvalues.  The  appropriate  boundary  conditions  are  E_  ~  r  near  the 
origin  and  E  =  0  at  r  =  a  (one  solution  of  Eq.  37  can  be  shown  to  have  such 
behavior  at  r  =  a).  One  begins  the  integration  at  r=0  and  adjusts  the 
eigenvalue  A  in  successive  integrations  so  as  to  assure  the  the  eigenfunction 
vanishes  at  r=a.  The  radial  dependences  of  the  first  few  m=2  modes  obtained 
in  this  manner  are  shown  in  Fig.  4.  One  can  see  from  Fig.  4  that  the 
turning  point  for  successively  higher  radial  trapped  Bernstein  modes  is  closer 
to  the  plasma  edge  and  would,  for  some  higher  radial  mode  number,  approach  the 
wall.  Under  these  circumstances,  the  evanescent  layer  becomes  thin  and  one 
can  no  longer  argue  that  the  conditions  at  the  edge  of  the  plasma  are 
unimportant.  Thus  this  approximate  approach  breaks  down  for  very  high  radial 
mode  numbers.  The  eigenvalues  obtained  are  plotted  in  Fig.  5  for  the  range 
0  <  p/a  <  .05.  For  very  small  values  of  this  parameter  (low  temperature)  the 


Bernstein  mode  spectrum  is  very  dense,  approaching  a  continuous  spectrum,  as 
found  earlier,  when  T  O.  For  a  given  m  >  0,  the  eigenvalues  all  have  A  < 
m-1,  and  for  this  profile  A  >  (m-l)/2. 


VIII.  EXTENSION  TO  MULTIPLE  ION  SPECIES.  Only  a  few  changes  are  required  to 
extend  the  previous  results  to  two  or  more  ion  species.  This  extension  is 
motivated  by  the  recent  experiments  of  Sarid  et  al^^’  First,  we  limit  the 
discussion  to  two  species  for  simplicity.  Second,  because  the  particles  have 
positive  charge,  they  gyrate  in  the  opposite  direction,  the  plasma  rotation  is 
in  the  opposite  direction,  and  fields  which  can  drive  cyclotron  resonance  have 
the  opposite  polarization.  In  the  discussion  below,  we  simply  ignore  these 
changes  in  sign  and  extend  the  previous  analysis  for  electrons  of  different 
mass,  and  make  appropriate  changes.  Thus  the  quantities  and  below  are 
positive.  As  before  we  characterize  the  density  profiles  of  the  two  species 
by 

noi(r)  =  n^^(O)  f^(r)  f^(0)  =  1  [40a] 

The  E  X  B  angular  velocity  of  each  species  is  the  same  so  we  define  a  common 
angular  velocity  function  g(r), 


where 


a;o(r) 


g(r) 


g(0)  =  1 


^oCO) 


2(jJ 


cl 


2(jj 


c2 


[41] 

[42] 


2  2  2 

If  we  let  ^^p2^^^c2  ^  ^pl^^^cl^  charge  fraction  of 

species  2,  then  Eq.  8a  becomes  (using  a  similar  definition  for  e^. 


When  f  (r)  and  f  (r)  are  specified  the  g(r)  is  determined.  For  example  if 
1  2 

f  (r)  =  1  -  a  (r/a)^  and  f-(r)  =  1  -  a  (r/a/)^  then 
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,  .  “i^i  “2^2  ,r^2 

B(r)  =  1 - 2 - 


[441 


Parabolic  profiles  serve  as  simple  illustrative  examples  of  diffusion-like 
non-rigid  rotor  profiles,  and  reduce  to  T=0  rigid  rotor  profiles  when  the 
a’s  are  set  to  zero. 

We  now  examine  cyclotron  resonance  when  the  frequency  is  near  the 

cyclotron  frequency  of  species  1.  The  dynamics  of  species  2  can  be  ignored, 

since  it  is  far  from  resonance.  However  its  effect  on  has  been  included 

in  g(r).  As  before,  we  introduce  a  normalized  frequency  shift,  ,  for  the 

mode  near  u  fn  u  normalized  to  the  central  rotation  frequency 
Cl , 


u  =  (jj  +  w  (0)  A 

Cl  ^  1 


[45] 


The  analysis  leading  up  to  Eq.  19  is  unchanged,  except  that  f(r)  is  replaced 
by  The  dielectric  function  for  circularly  polarized  fields  becomes 


K(r)  =  1  - 


'1 


A^-(m-l)g(r)  -.5rdg(r)/dr 


and  when  we  use  Eq.  43  and  solve  for  1/K  -1  we  obtain 

1 


K(r) 


-  1  = 


Ar(m-l)g+S^2  ■ 


[46] 


with  a  similar  expression  (with  subscripts  1  and  2  interchanged)  when  u  ^^2* 
For  the  plasma  extending  to  r=a  and  the  wall  at  r=b  we  get 


m  b^^ 


■ 


-  2m-lj 
2mr  dr 


[47] 


This  gives  the  behavior  when  the  frequency  is  in  the  vicinity  of  and  a 


similar  expression  (with  the  subscripts  1  and  2  interchanged)  gives  the 
behavior  when  the  frequency  is  in  the  vicinity  of  u  .  When  e  =0,  no  second 
species,  and  m=l  we  recover  Eq.  33,  our  previous  single  species  result.  As  in 
the  single  species  case,  the  two-species  integrals  such  as  Eq.  47  can  be 
evaluated  in  the  case  of  a  parabolic  density  profile,  and  this  evaluation 
leads  to  logarithmic  functions  as  given  in  Eq.  35.  A.  Plasma  Density 
Independent  of  Radius.  For  rigid  rotor,  or  "top  hat"  profiles, 
f  (r)=f  (r)=g(r)=l  for  r^a  so  that 

L  Z 


R„(b)  =  - 
m 


_  1  / a  N  2m 

A  -(m-l)+e  b  ’  ^  ^  ^cl 

X  z 

e 


R  (b)  =  —  "T  7  -  VTT 

m  X  -{m-l)+€  b 

A  1 


2  2m 


(jJ  pa  tJ 


c2 


The  resonant  modes  are  obtained  by  setting  equal  to  unity  and  this  gives 


\  -  S 

\  =  (m-1)  - 


2m 


a _ 

b2m 

2m 


a _ 

b2m 


a;  w  w  ,  [48a] 

cl 

OJ  pa  u  ^  [48b] 


Upon  comparison  with  Eq.  32  we  see  that  the  terms  involving  e  describe  the 

changes  in  resonant  frequencies  arising  from  the  presence  of  a  another 

species.  For  the  center  of  mass  modes  (m  =  1), 

2  2 


A 


1 


[49a] 


A 


2 


1  . 


[49b] 


where  we  have  used  the  fact  that  e  +6  =1 .  The  second  term,  proportional  to 

X  z 

the  density  of  the  other  species,  represents  an  additional  downshift  caused  by 
the  presence  of  the  other  species. 


B.  Plasma  Density  Dependent  Upon  Radius.  We  now  consider  density  profile 
effects  and,  for  simplicity,  limit  our  discussion  to  the  parabolic  profiles 
used  earlier.  When  we  examine  Eq.  47,  and  the  similar  expression  with 
subscripts  1  and  2  interchanged,  we  note  that  it  is  possible  for  the 
denominator  of  the  integral  to  vanish  at  some  radius  in  certain  frequency 
bands.  In  the  cold  plasma  theory  this  signifies  absorption  bands  and  in  the 
warm  plasma  theory  these  are  the  bands  in  which  the  Bernstein  modes  may  be 
found.  To  determine  these  frequency  bands  consider  the  vanishing  of  the 
denominator  with  the  assumed  profiles.  Note  that  bands  associated  with  each 
of  the  species  is  expected.  The  denominator  for  the  species  1  resonance  has 
the  form 

A  .(„-l)(l-|4)  «,(!-«  4* 

a  a 

where,  for  simplicity,  we  assume  that  a  =  a  =  a,  i.e.  that  both  species  have 
the  same  density  profile.  Vanishing  of  this  denominator  at  r=0  and  r=a 
serves  to  determine  the  limits  of  the  absorption  bands.  These  limits  are, 
for  species  1, 

A,  =  (m-l)-£,  [50a] 

1  2 

Aj  =  (m-l)(l-|)- 

with  similar  expressions  for  the  limits  associated  with  the  second  species. 
For  the  center  of  mass  modes  (m=l),  the  frequencies  of  these  two  bands  can  be 
described  as 

-fj  ^  A^  ^  -e^Cl-a)  [51a] 

^  Aj  ^  -fj(l-a),  [51b] 

i.e.  both  bands  lie  below  their  respective  cyclotron  frequencies.  From  this 
result,  we  see  that,  unless  the  density  is  independent  of  radius  (a=0. 


top-hat  distributions),  there  are  now  absorpt ion  bands  associated  with  the 

in=:l  cyclotron  modes,  whereas  for  the  single  species  case,  there  was  only  a 

single  discrete  resonance.  We  speculate  that  the  inclusion  of  temperature 

will  lead  to  Bernstein  modes  in  these  bands. 

IX.  SUMMARY  AND  CONCLUSIONS. 

We  have  studied  the  k^=0  cyclotron  modes  of  a  cylindrical  nonneutral 

2  2 

plasma  of  low  density  (w  /w  «  1).  Particular  attention  has  been  given  to  a 

P  c 

single  species  (electrons)  which  is  in  nonuniform  rotation,  i.e.  where  the 
density  and  angular  velocity  are  functions  of  radius.  We  have  illustrated  the 
general  features  of  the  behavior  near  cyclotron  resonance  by  using  a  simple 
parabolic  profile.  The  qualitative  behavior  for  other  monotonical ly 
decreasing  profiles  is  similar  but  differs  in  details.  A  We  have  obtained  the 
admittance  function,  Y(a;),  which  characterizes  the  current  flowing  to  a 
segment  of  the  wall  when  an  oscillating  potential  whose  frequency  is  near  the 
cyclotron  frequency  is  applied  to  that  segment.  An  absorption  band  is  found 
for  each  m  value,  for  the  range  of  frequencies  for  which  a  local  upper  hybrid 
resonance  occurs  at  some  radius  in  the  nonuniform  plasma.  These  absorption 
bands  correspond  the  the  frequency  bands  in  which  the  plasma  has  a  continuum 
of  singular  eigenfunctions,  although  that  approach  is  not  persued  here.  The 
plasma  responds  only  to  the  circularly  polarized  component  of  the  electric 
field  which  has  the  same  sense  as  the  gyrating  charged  particles  and,  because 
of  plasma  rotation,  the  each  absorption  band  is  upshifted  by  the  Doppler 
effect.  m=l  is  found  to  be  a  special  case  with  only  a  single  discrete  mode 
and  no  continuum,  a  property  which  it  shares  with  the  low  frequency  m=l 
diocotron  mode.  Both  are  center  of  mass  modes.  Presumably  in  both  cases  an 
m=l  continuum  of  singular  eigenmodes  exists  but  they  have  no  electric  field 


exterior  to  the  plasma,  and  therefore  not  observable,  or  excitable,  with 
electrodes  exterior  to  the  plasma.  The  m=l  cyclotron  mode  was  shown  to  have  a 
frequency  which  lies  below  the  cyclotron  frequency  by  and  amount  equal  to  the 
low  frequency  diocotron  mode  frequency  <w  >/2w  ,  where  the  average  is  taken 
over  the  cross  section  of  the  cylinder. 

When  finite  Larmor  radius  effects,  owing  to  plasma  temperature,  are  taken 
into  account  in  an  approximate  fashion,  Bernstein  modes  are  found  in  the 
frequency  bands  for  which  the  cold  plasma  theory  predicts  an  absorption 
continuum  associated  with  the  continuus  spectrum  of  singular  eigenmodes. 
Because  of  the  radial  density  and  angular  profiles  the  Bernstein  modes  are 
trapped  in  the  central  core  of  the  plasma.  They  are  traveling  waves  in  the 
azimuthal  direction  and  standing  waves  in  the  radial  direction.  The  relative 
frequency  spacing  between  modes  with  different  numbers  of  radial  nodes 
is  proportional  to  the  ratio  of  Larmor  radius  to  gradient  scale  length,  a 
small  quantity.  The  Bernstein  mode  spacing  can  be  used  to  estimate  the  plasma 
temperature.  Finally,  we  have  outlined  an  extension  of  the  cold  plasma 
results  to  plasmas  with  multiple  ion  species. 

A  more  complete  theory  of  the  trapped  Bernstein  modes  in  a  cylindrical 
nonuniform  nonneutral  plasma,  such  as  a  solution  of  the  Vlasov  equation  to 
first  order  in  the  temperature,  is  desirable.  This  approach  has  been  used  in 
Ref.  6  for  a  neutral  plasma  and  in  Ref.  8  for  a  planar  nonneutral  plasma.  In 
order  to  carry  out  such  an  analysis  one  would  also  have  to  find  a  suitable 
steady  state  with  radial  profiles  of  the  density  and  angular  velocity  which 
includes  the  diamagnetic  drifts  as  well  as  the  E  x  B  drifts. 
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Schematic  diagram  of  the  plasma,  surrounding  electrodes,  and 

driving  circuitry  (end  view),  with  definition  of  Y(w) ,  the  plasma 

admittance  function  with  u>  =  u  +  w  (0)A. 

c  0 

Illustrative  density  and  angular  velocity  profiles  of  a  rotating 
nonneutral  plasma,  normalized  to  their  central  values. 
f(x)  =  np(x)/n^(0),  g(x)  =  w^(x)/u;^(0),  x  =  r/a. 

Admittance  versus  frequency  for  the  m  =  1,2, 3, 4  modes  and  the 

profiles  shown  in  Fig.  2.  Plotted  are  the  real  and  imaginary 

parts  of  the  plasma  contribution,  2R  (b)/[R  (b)-l].  The  ratio  of 

m  m 

wall  radius,  b,  to  plasma  radius,  a  is  1.25.  Plots  for  m  =  3  and 
4  have  been  expanded  by  factors  of  5  and  10  with  respect  to  m  =  2 
for  clarity  in  presentation. 

Radial  Wave  Functions  for  m=2,  p/a=0 .02,  with  their  normalized 
frequencies.  The  first  five  radial  modes  are  shown. 

Theoretical  spacing  of  the  trapped  m=2  Bernstein  modes  versus 
Larmor  radius,  p/a,  normalized  to  the  central  rotation  frequency 


FIG.  1 .  Schematic  diagram  of  the  plasma,  surrounding  electrodes,  and  driving  circuitry  (end 
view),  with  definition  of  Y((0),  the  plasma  admittance  function  with  0)  =  (0^.  +  (Oo(O)X  . 


FIG.  2.  Illustrative  density  and  angular  velocity  profiles  of  a  rotating  nonneutral  plasma, 
normalized  to  their  central  values:  f(x)  =  no(x)/no(0),  g(x)  =  (q^(x)/(S)J(0),  x  =  r/a. 
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FIG.  3.  Admittance  versus  frequency  for  the  m=l,2,3,4  modes  and  the  profiles  shown  in  Fig.  2. 
Plotted  are  the  real  and  imaginary  parts  of  the  plasma  contribution,  2R^(b)/[R^(b)-l].  The  ratio 
of  wall  radius  to  plasma  radius,  b/a,  is  1 .25.  Plots  for  m=3  and  4  have  been  expanded  by  factors 
of  5  and  10  with  respect  to  the  m=2  for  clarity  in  presentation 


Fig.  4.  Radial  wave  functions  for  m  =  2,  p/a  =  0.02,  with  their  normalized  frequencies. 
The  first  five  modes  are  shown. 


Fig.  5.  Theoretical  spacing  of  the  trapped  m  =  2  Bernstein  modes  versus  Larmor  radius,  p/a, 
normalized  to  the  central  rotation  frequency  ^^l  2f^  for  the  profiles  of  Fig.  2. 
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8R  19  Negative  Energy  and  Dapetaioa  of  the  Dioeotioa  Reaonanwa*. 
Sateesh  Pillai  and  Roy  W.  Gould,  (Mifonua  •/  Teeknoton.  A 

cylindrical  column  of  pure  electron  plasma  exhibits  characteristic  quasi¬ 
mode  “diocotron"  resonances  around  the  angular  rotation  frequencies. 
These  are  resonances  caused  by  the  fluid  nature  of  the  qrstem  and  are 
identified  with  the  eigen  number  *m”  corresponding  to  iU  circular 
symmetry.  It  has  been  experimentally  demonstrated  that  the  undamped 
m=l  diocotron  is  a  negative  energy  mode*.  The  m=2  quasi-mode  is 
normally  damped  by  a  collisionless  process^ .  We  present  first 
experimental  evidence  showing  that  the  addition  of  external  resistors  to 
the  wall  segments  decrease  the  damping  rate  for  the  m=2  mode, 
demonstrating  that  the  m=2  diocotron  is  also  a  negative  energy  mode. 
Based  on  the  linear  theory  the  higher  diocotron  modes  are  also  negative 
energy  modes. 

The  dispersion  characteristics  of  the  diocotron  modes  for  a 
particular  profile  will  also  be  presented. 
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8R20  Numerical  Collisionless  Damping  of  the  Diocotron 
Resonance*.  David  A.  Bachman  and  Roy  W.  Gould,  California 
InslUule  of  Technology.  A  ^lindrical  pure  electron  plasma  features 
diocotron  resonances  near  die  angular  rotation  frequen^.  These 
resonances  are  characterized  1^  dieir  circular  symmetry. 
Experimentally,  the  m=2  diocotron  mode  is  found  to  dec^  by  a 
collisionless  process.  Additionally,  a  solution  of  die  linearized 
fluid  equations  for  representative  density  and  angular  profiles  also 
gives  collisionless  deci^. 

A  two  dimensional  nonlinear  fluid  code  of  a  pive  electron 
plasma  has  been  used  to  study  die  plasma  response  to  an  applied 
pulse.  The  fluid  code  also  exhibits  the  collisionless  damping  of  die 
m-2  diocotron  mode,  and  die  damping  rate  of  diis  mode  found 
from  the  nonlinear  fluid  code  matches  die  damping  rate  found 
from  the  linearized  numerical  calculation. 
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